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ABSTRACT 

This  report  was  prepared  to  describe  the  results  of 
an  experimental  study  of  free  turbulent  mixing  in  axial 
pressure  gradients.   Both  the  pressure  field  and  the  velocity 
field  were  extracted  using  pitot-static  probes  and  both  a 
single  normal  wire  and  an  X-array  hot-wire  anemometers. 
The  resulting  hot-wire  data  shown  in  Figures  1-80  present 
an  excellent  visual  representation  of  the  physical  aspects 
of  turbulent  flow. 

An  analysis  of  the  results  of  this  experiment  show 
conclusively  that  the  Reynolds   Stresses  in  the  flow  are 
related  by  some  function  to  the  turbulent  kinetic  energy; 
however,  no  specific  function  was  found  that  would  produce 
the  relationship  observed  in  this  study. 
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NOMENCLATURE 

A  Boussinesq  apparent  viscosity 

A  X-Array  hot-wire  lead  (u+v) 

B  X-Array  hot-wire  lead  (u-v) 

C  Constant  in  hot-wire  calculations 

C,  Constant  in  apparent  viscosity  relation 

C  Local  turbulent  shear  coefficient 

T 

C  Correlation  coefficient 

d  Wake  generator  diameter 

e  Linear  output  from  hot-wire 

H  Shape  factor 

l,j  Indicates  i  or  j  direction  (1=1,2,...   j-1,2,...) 

k  Turbulent  kinetic  energy 

I  Prandtl  mixing  length 

p  Pressure 

q  Dynamic  pressure 

r  General  flow  parameter 

S(x)  Apparent  viscosity  constant 

T  Characteristic  time  / 

t  Time 

x,y,z  Spatial  Coordinates 

x  -  Axial 

y  -  Lateral 

u,v,w  Corresponding  velocities 

v*  Characteristic  velocity 


B  Constant  in  hot-wire  calculations 

6  Displacement  thickness 

6-,  Energy  thickness 

e  Mixing  length  parameter  =  l^sP 

0  Momentum  thickness 
u  Laminar  viscosity 
u.  Apparent  viscosity 

v  Laminar  kinetic  viscosity 

v.  Apparent  kinetic  viscosity 

p  Density 

t  Reynold's  Stress 

V      Differential  operator  DEL 
b 
f   dy   Integral  from  a  to  b  with  respect  to  y 

a 

/ 

8u 

7r—     Partial  derivative  with  respect  to  r 

■=■—     Ordinary  derivative  with  respect  to  r 

Indicates  absolute  value 

Subscripts 

A      Parameter  of  lead  A  on  hot-wire 

B      Parameter  of  lead  B  on  hot-wire 

e      Local  external  flow  conditions 

x,y    Indicates  partial  differentiation  with  respect  to  x,y 

t      Shear  related  quantity 

Superscripts 

Indicates  a  time  averaged  quantity 

1  Designates  a  fluctuating  quantity 


I.   INTRODUCTION 

In  most  real  fluid  flows  the  transition  from  laminar  to 
turbulent  flow  carries  the  fluid  mechanist  from  a  situation 
where  he  can  safely  predict  the  values  of  the  flow-field 
parameters  from  the  solution  of  the  Navier-Stokes  Equations 
to  one  where  predictions  are  obtainable  only  on  the  average, 
at  best.   Yet  most  practical  applications  require  at  least 
some  consideration  of  turbulence  and  how  it  can  be  dealt 
with.   In  the  past  the  engineer  has  assumed  various  plausible 
hypotheses  in  order  to  introduce  a  measure  of  order  and 
predictability  to  turbulent  flow.   It  is  the  intention  of 
this  study,  then,  to  attempt  to  experimentally  verify  certain 
of  these  hypotheses.   Before  describing  the  currently 
accepted  assumptions  regarding  turbulent  flow,  a  general 
review  of  their  history  will  be  presented  to  provide  the 
necessary  background  information. 

A  turbulent  flow  is  defined  as  one  in  which  an  irregular 
fluctuation  or  mixing  motion  is  superimposed  on  the  princi- 
pal fluid  motion.   Its  onset,  is  accompanied  by  a  noticeable 
change  in  the  law  of  resistance.   In  the  case  of  laminar 
pipe  flow,  the  longitudinal  pressure  gradient  which  main- 
tains the  motion  is  proportional  to  the  first  power  of  the 
velocity.   By  contrast,  in  turbulent  flow  this  pressure 
gradient  becomes  nearly  proportional  to  the  square  of  the 
mean  flow  velocity.   In  order  to  calculate  flow  field  values 


in  a  turbulent  field,  then,  some  logical  hypothesis  must 
be  made  to  account  for  this  change  in  the  apparent  viscosity 
of  the  flow.   First,  it  is  generally  assumed  that  the  fluid 
motion  can  be  separated  into  a  mean  and  a  fluctuating  part 
which  are  linearly  related  to  the  total  fluid  motion  in 
such  a  way  that  a  flow  parameter  (r)  can  be  defined 

r  =  r  +  r'  (1-1) 

These  definitions  can  be  substituted  into  the  equations  of 
motion  to  produce  a  new  set  of  equations  involving  mean  and 
fluctuating  terms  and  their  cross  products.   Assuming  that 
the  time  average  of  the  fluctuations  is  zero,  the  equations 
may  be  averaged  to  get  the  following  momentum  equation  for 
a  two-dimensional  mean  flow  (see  Appendix  A  for  the  complete 
development ) : 


<■> ~ dll  r\  2  — 

— dU  .  -  dU  e         /sdU       fr     ot\ 

"^  +  v  w  =  ue  ar- +  vt  (x)  7-2       (A~22) 

ay 


where  the  apparent  kinematic  viscosity  (v, (x))  is  a  function 
of  the  turbulent  flow  field,  and  is  calculated  by  using  some 
logical  empirical  hypothesis  for  approximating  its  value. 
Since  only  the  mean  values  of  the  flow  parameters  are 
known,  some  means  of  relating  these  values  to  the  Reynolds 
stresses  produced  by  the  fluctuating  part  of  the  flow  is 
required.   The  first  hypothesis  for  relating  these  quantities 


was  proposed  by  J.  Boussinesq, 


T   =  -pUMTT  =  A   |u  (1-2) 

o  t  3y 


where  A  is  an  apparent  or  eddy  viscosity  that  is  not  a 
fluid  property  but  depends  on  the  structure  of  the  flow 
field. 


At  =  p  vt(x) 


This  apparent  viscosity  is  generally  assumed  to  be  propor- 
tional to  some  characteristic  velocity  (v1)  times  the  mean 
velocity  gradient: 


AT  -  v*  |f  (1-3) 


In  1925  Prandtl  made  an  important  advance  in  turbulent 
theory  by  introducing  the  concept  of  mixing  length  and  the 
relation, 

x   =  p  £2  jdU|  d^  (I_4) 

t    p     'dy1  dy 

v,  (x)  =  I        -3— 
t         '  dy  ' 


where 


£  =  mixing  length 


du 


^  \tt~ I  =  characteristic  velocity 


This  relation  has  produced  relatively  accurate  values  for 
some  limited  cases.   However,  this  definition  of  the  kine- 
matic viscosity  implies  that  it  should  vanish  where  the 
velocity  gradient  is  zero,  which  is  contradictory  to  experi- 
mental evidence.   Several  subsequent  hypotheses  were  proposed 
by  Taylor,  Von  Karman,  and  Prandtl  himself,  but  they  were 
all  limited  in  application  to  free  turbulent  flow. 

Later,  in  an  effort  to  develop  a  universally  valid 
model  for  the  function  y , ,  Prandtl  introduced  an  energy- 
based  hypothesis  which  was  subsequently  improved  on  by 
Emmons  (195*0  >  and  Glushko  (1965), 


Wt  «  p  k1/2  I  (1-5) 


k  s  (u2  +  v"2)/2 


Here  k  was  found  to  be  a  better  characteristic  velocity  than 
in  the  earlier  mixing  length  hypothesis,  but  the  length 
scales  were  difficult  to  ascertain.   More  modern  theories, 
referred  to  as  two-equation  viscosity  models,  use  sophisti- 
cated analytical  methods  to  determine  this  length  scale. 
The  apparent  viscosity  is  first  assumed  to  be  related  to  the 
kinetic  energy.   Then  some  hypothesis  of  the  form 


e  =  km  ln 
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is  made  to  make  it  possible  to  calculate  the  characteristic 
length.   One  of  the  most  accepted  of  these  relations  was 
introduced  by  Chou  in  19^5: 


e  =  k3/2/l  (1-6) 


2 
yt  =  C1  p  k  /e 


Thus 


yt  =  C1  P  k1/2  £  (1-7) 


It  was  this  type  of  assumption  which  initiated  the  pro- 
ject reported  here.   Using  both  pitot-static  and  hot-wire 
anemometer  sensors,  an  attempt  was  made  to  validate  the 
hypothesis  that  shear  was  predominantly  a  function  of  the 
turbulent  kinetic  energy. 
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II.   EQUIPMENT 

A.   WIND  TUNNEL 

1 .  General  Description 

The  tunnel  consists  essentially  of  three  major 
components:  the  inlet,  the  test  section,  and  the  blower  as 
shown  in  Figures  81-83.   The  inlet  was  designed  to  contract 
the  slow  moving  inlet  air  to  a  high-velocity  parallel  flow 
at  the  entrance  of  the  test  section.   At  this  point,  a 
plexiglas  cylinder  was  introduced  into  the  flow  to  generate 
a  wake  velocity  defect  which  was  then  subjected  to  a  constant 
adverse  pressure  gradient.   After  leaving  the  test  section, 
the  flow  was  maintained  at  a  constant  area  through  a  length 
of  tubing,  then  diffused  and  directed  into  the  centrifugal 
blower  and  exhausted  vertically. 

2.  The  Inlet 

The  inlet  section  of  the  tunnel  was  composed  of  the 
inlet  bellmouth,  the  plenum  chamber,  and  the  nozzle.   Incoming 
air  was  first  contracted  into  a  two-foot  square  plenum  sec- 
tion.  This  contraction  was  accomplished  by  the  inlet  bell- 
mouth  constructed  from  eight-inch  diameter  aluminum  pipe 
halves.   The  corner  seams  were  filled  with  a  sanded  epoxy 
fillet  to  effect  a  smooth  transition  of  the  air  into  parallel 
flow.   The  plenum  section  consisted  of  four  two-foot  square 
by  one-foot  long  steel  sections  bolted  together.   Cork  gas- 
kets were  inserted  between  the  sections  to  prevent  the  influx 
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of  air  at  the  joints  and  to  allow  for  undisturbed  air  flow 
over  these  junctions.   A  four-inch  thick  section  of  quarter- 
inch  honeycomb  was  installed  to  straighten  the  air  flow  and 
prevent  swirling.   Next,  two  copper  screens  were  inserted 
to  reduce  the  ambient  turbulence  level  in  the  flow.   Finally, 
the  air  passed  from  the  plenum  chamber  to  the  test  section 
through  a  two-foot  long  steel  nozzle  with  an  exit  area  of 
forty-nine  square  inches  (seven  inches  square),  and  produced 
a  contraction  ratio  of  11.76.   Static  pressure  ports  were 
located  on  the  right  side  of  the  nozzle  at  the  inlet  and 
exit  for  determination  of  the  exit  dynamic  pressure. 
3.   The  Test  Section 

The  test  section  consisted  of  two  parallel,  vertical 
side  walls  separated  by  flexible  and  adjustable  top  and 
bottom  walls .   The  flexible  walls  were  constructed  from 
one-eighth  inch  Plexiglas  fastened  to  the  nozzle  with  piano 
hinges.   As  described  in  Reference  1,  these  walls  could  be 
varied  to  produce  a  variety  of  different  pressure  gradients, 
but  for  this  experiment  they  were  set  to  create  a  constant 
adverse  pressure  gradient  followed  by  a  constant  favorable 
pressure  gradient. 

One  vertical  wall  was  provided  with  thirteen  slots 
for  inserting  data  sensors.   These  slots  were  nineteen 
inches  high,  three-eighths  of  an  inch  wide,  and  spaced  five 
inches  apart.   A  filler  plug  for  each  slot  was  carefully 
hand  milled,  and  polished  to  prevent  disturbing  the  flow 
along  the  wall  when  the  slots  were  not  in  use.   A  thin 
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aluminum  strip  was  wedged  over  each  filler  plug  to  ensure  a 
tight  fit.   When  these  plugs  were  removed  to  insert  probes, 
a  sliding  teflon  seal  was  used  to  prevent  an  influx  of  air 
through  the  slot  (see  Figure  8*0. 

The  other  vertical  wall  of  the  test  section  had  three 
rectangular  access  holes  along  the  axial  centerline  for 
adjustments  and  measurements.   In  addition,  as  one  can  see 
in  Pig.  85,  thirty-three  static  pressure  ports  were  inserted 
along  the  axial  centerline  at  one-and-one-half  inch  inter- 
vals.  A  manometer  bank,  shown  in  Fig.  86,  connected  to  these 
static  pressure  ports  provided  a  visual  indication  of  the 
uniformity  of  the  pressure  distribution,  as  well  as  indica- 
ting the  presence  of  any  regions  of  separation.   As  an 
additional  indication  of  possible  flow  separation,  several 
rows  of  tufts  were  attached  to  the  walls  in  critical  areas. 

Having  produced  the  desired  environment  for  the 
type  of  air  flow  that  this  study  was  initiated  to  investi- 
gate, a  turbulent  wake  velocity  defect  was  introduced  into 
the  tunnel  by  placing  a  one-inch  diameter  cylinder  in  the 
inlet  of  the  test  section.   In  order  to  ensure  a  turbulent 
boundary  layer  at  separation  from  the  cylinder,  a  thin  strip 
of  abrasive,  shown  in  Fig.  87,  was  attached  to  the  stagnation 
line  of  the  cylinder. 
iJ.   The  Blower 

After  exiting  the  test  section,  the  air  passed 
through  two  seven-inch  square,  four-foot  long  steel  sections 
into  a  diffusor  and  then  into  the  blower.   A  centrifugal 
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blower  operating  at  three  thousand  revolutions  per  minute 
was  utilized  to  pull  air  through  the  tunnel.   The  air  was 
then  exhausted  vertically  through  k    sections  of  steel  pipe 
which  contained  two  dampers  for  the  control  of  the  total 
mass  flow  through  the  tunnel. 

B.   DATA  SENSORS 

1.  General  Description 

The  data  for  this  study  were  taken  with  the  three 
different  types  of  probes  shown  in  Fig.  88:   A  pitot  static 
probe,  a  single  normal  hot  wire  probe,  and  an  x-wire  hot 
wire  probe.   A  detachable  traversing  mechanism  was  used  to 
position  the  probes  vertically,  and  the  values  were  recorded 
by  hand  or  on  an  x-y  plotter. 

2 .  The  Traverse 

This  apparatus  was  equipped  with  two  constant  speed 
motors  to  drive  the  mechanism  across  the  nineteen-inch 
vertical  span  of  the  tunnel.   The  probe  was  clamped  to  a 
movable  support  on  the  traverse,  and  the  vertical  position 
measured  by  indexing  a  reference  mark  on  this  mechanism  to  a 
scale  fixed  to  the  frame  of  the  apparatus.   This  scale  was 
readable  to  one  one-hundreth  of  an  inch.   A  rack-and-pinion 
gearing  system  drove  a  potentiometer  that  was  also  attached 
to  the  movable  support  to  produce  an  electrical  position 
signal  for  the  abscissa  of  the  x-y  plotter. 
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3.  The  Pitot-Static  Sensor 

A  Prandtl-type  pitot-static  pressure  probe  Model 
#   PBA-12-F-11-KL,  built  by  United  Sensors  and  measuring 
8l  hundreths  of  an  inch  from  stem  to  tip,  was  used  to  obtain 
dynamic  and  local  static  pressures  at  one-tenth  of  an  inch 
intervals  through  the  wake .   These  pressures  were  read  on 
two  variable  angle,  water-filled  U-tube  manometers  referenced 
to  atmospheric  pressure.   This  reference  was  chosen  to  re- 
duce the  fluctuations  in  the  water  levels.   As  described  in 
Reference  1,  taking  P-total  minus  P-static  directly  creates 
very  large  fluctuations  that  make  accurate  data  taking  diffi- 
cult.  By  separating  the  two  readings,  the  fluctuations  in 
the  manometer  were  significantly  reduced,  but  not 
eliminated. 

A  third  manometer  was  attached  to  the  plenum  chamber 
of  the  tunnel,  and  referenced  to  the  last  static  pressure 
port  on  the  test  section.   As  one  can  see  in  Fig.  89,  all  of 
these  manometers  were  equipped  with  portable  plastic  sights 
which  restricted  one's  line  of  sight  to  a  plane  perpendicu- 
lar to  the  manometer,  thus  precluding  possible  parallax  and 
assuring  consistent  readings. 

4.  The  Normal  Wire 

This  sensor,  shown  in  Fig.  90,  was  a  single-wire 
anemometer  probe,  Model  #   1150  built  by  Thermo-Sy stems  Inc. 
(TSI),  measuring  4  inches  from  stem  to  tip.   Again  the  probe 
was  fixed  to  the  traverse  for  data  collection.   The  lead 
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from  the  sensor  was  connected  to  a  105^A-30  anemometer 
module  with  a  1051-2  monitor  and  power  supply,  also  built 
by  Thermo-Systems  Inc.  (TSI)  (See  Fig.  91  &  92),  which 
converted  the  voltage  changes  to  readable  data.   Next  these 
data  were  fed  into  a  TSI  1015C  correlator  that  was  used  to 
separate  the  alternating  and  D.C.  signals  (see  Appendix  B). 
The  output  of  the  correlator  was  connected  to  an  analog 
data  processing  computer  which  indicated  the  mean  square 
average  output  of  the  correlator.   Finally,  the  output  of 
this  apparatus  was  connected  to  the  ordinate  of  the  X-Y 
plotter,  recalling  that  a  potentiometer  on  the  traverse  was 
connected  to  the  abscissa. 

5.   The  X-Wire  Anemometer  Probe 

The  last  sensor,  a  Model  #  1055  built  by  TSI,  con- 
sisted of  two  hot  wires  oriented  at  45°  to  the  axis  of  the 
probe  such  that  the  output  signals  of  the  wires  corresponded 
to  u  +  v  and  u  -  v  respectively  (See  Appendix  B) .   These 
data  were  processed  in  the  same  manner  as  the  normal  wire 
data  except  for  the  use  of  several  additional  algebraic 
functions  of  the  1015C  correlator. 
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III.   EXPERIMENTAL  PROCEDURE 

A.   GENERAL  DESCRIPTION 

The  previous  study  described  in  Reference  1  had  con- 
figured the  test  section  to  produce  a  constant,  adverse- 
favorable  pressure  gradient.   Since  this  was  the  environment 
required  for  this  study,  no  changes  were  made  in  the  test 
section  except  to  ensure  a  tight  seal  between  the  vertical 
and  the  movable  walls  by  the  addition  of  strips  of  clay 
along  the  lines  of  contact. 

The  actual  collection  of  data  was  carried  out  in  three 
distinct  steps.   First,  the  pitot-static  system  was  used  to 
collect  the  total  and  static  pressure  data  that  were  intro- 
duced into  the  computer  program,  described  in  Appendix  E, 
and  used  to  calculate  the  flow  parameters.   Next,  a  normal- 
wire  anemometer  was  inserted  in  the  test  section  to  gather 
information  on  both  the  mean  and  fluctuating  axial  velocities 
These  hot-wire  values  were  converted  to  velocities  by  equa- 
ting the  maximum  voltage  at  each  station  to  the  maximum  value 
of  the  velocity  calculated  from  the  pitot-static  data. 
Finally,  an  X-array  hot-wire  anemometer  probe  was  used  to 
get  both  mean  and  fluctuating  data  for  velocities  in  both 
the  axial  and  the  lateral  directions,  as  well  as  Reynolds 
stresses.   Throughout  the  experiment  a  constant  wind  tunnel 
speed  was  maintained  by  ensuring  that  the  total  pressure  of 
the  Plenum  chamber,  minus  the  static  pressure  from  the  test 
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section,  was  no  less  than  9.90  nor  greater  than  10.00 
inches  of  water  on  a  thirty-degree  inclined  manometer. 

B.   PITOT-STATIC  DATA 

In  this  phase  of  the  experiment  all  data  were  recorded 
(to  an  estimated  accuracy  of  one  one-hundreth  of  an  inch) 
from  water-filled  U-Tube  manometers  connected  to  the  pitot- 
static  probe.   At  each  of  eleven  stations  downstream  from 
the  cylinder  the  static  and  total  pressures  were  recorded 
at  intervals  through  the  wake  of  one-tenth  of  an  inch.   At 
each  station  the  probe  was  moved  near  the  upper  limit  of 
the  traverse  and  situated  in  the  free  stream,  then  the  tra- 
verse was  marched  across  the  flow,  recording  the  values  of 
total  and  static  pressure  at  each  increment.   The  initial 
determination  of  the  freestream  was  accomplished  by  observing 
that  the  manometer  becomes  almost  totally  insensitive  to 
changes  in  the  vertical  direction  in  the  undisturbed  flow. 
A  change  in  total  pressure  of  less  than  three  one-hundredths 
of  an  inch  for  three  consecutive  positions  was  taken  to  be 
the  edge  of  the  wake  defect. 

As  the  traverse  approached  the  center  of  the  wake  the 
intensity  of  the  turbulence  became  so  large  that  accurate 
readings  were  difficult  to  obtain  due  to  fluctuations  in  the 
water  level.   However,  by  observing  the  water  levels  for  a 
long  period  of  time  one  could  determine  an  average  height 
that  is  estimated  to  be  accurate  to  three  one-hundredths  of 
one  inch.   The  same  criterion  was  used  at  the  bottom  of  the 
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flow  channel,  as  at  the  top,  to  determine  when  the 
freestream  had  been  reached. 

C.   AXIAL  PLOW  ANEMOMETER  DATA 

The  probe  used  in  this  phase  was  a  single  normal  hot- 
wire anemometer  designed  to  produce  voltage  changes  corres- 
ponding to  velocity  fluctuations  in  the  flow.   When  the 
tungsten  wire  was  exposed  to  the  air  in  the  wind  tunnel  it 
was  cooled  by  the  flow  over  it,  thus  requiring  additional 
voltage  to  maintain  a  constant  temperature  in  the  wire .   It 
is  apparent,  then,  that  the  velocity  in  the  tunnel  is  pro- 
portional to  the  voltage  required  to  maintain  a  constant 
temperature  in  the  hot-wire.   The  probe  is  placed  in  a 
bridge  circuit  in  the  TSI  module,  and  its  output  linearized 
to  produce  a  linear  function  of  velocity  of  the  form: 

e  =  c  (u  +  u!)  (III-D 


where, 


and 


u  =  u  +  u' 

u  =  mean  flow  velocity 

u*  =  fluctuating  velocity 


T 
/   u1  dT  e  0  (III-2) 

0 


Therefore  if  the  DC  output  on  the  correlator  is  selected, 
the  output  would  be 
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e  =  c (u  +  u' ) 

However,  since  this  signal  was  a  strong  function  of  time, 
the  data  were  averaged  using  a  low  pass  filter  with  a  time 
constant  (T)  equal  to  0.6  sees.   The  final  output  of  the 
averager  was  simply 

e  =  cu  (III-3) 

and  this  quantity  was  plotted  versus  the  traverse  position. 

The  second  type  of  data  for  this  phase  was  taken  with 
the  selector  switch  of  the  correlator  on  AC,  thus  producing 
an  output  corresponding  to  the  fluctuating  portion  of  the 
flow: 

e'  =  cu'  (III-H) 

This  signal  was  also  a  strong  function  of  time  but,  because 
of  the  nature  of  the  term  u' ,  as  expressed  in  (III-2),  the 
signal  could  not  be  directly  averaged.   Therefore,  in  order 
to  provide  a  finite  measure,  u'  was  first  squared  and  then 
averaged  and  plotted  versus  the  traverse  position,  since 


i  /   (u')2  dx  t   0 
1  0 
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Once  a  plot  of  e  and  e'  had  been  obtained,  the  velocity 
could  be  calculated  by  solving  equations  (III-3)  and  (III-4) 
to  get 

u  =  1/C  e 


u'2  =  \   e'2 
C 


where  from  equation  (B-5)  in  Appendix  B 

.  e~~ 
C  E-i-  (HI-5) 

u 
e 

In  order  to  produce  accurate  values  for  u  and  u'  in  the 
laboratory  environment  the  equipment  had  to  be  carefully 
adjusted  for  each  individual  flow.   Specific  instructions 
for  this  procedure  can  be  found  in  the  TSI  operational 
manual  but,  in  general,  the  first  step  was  to  adjust  the 
control  resistance  to  balance  the  bridge  circuit..  This 
adjustment  was  accomplished  by  setting  the  control  resistor 
to  a  value  that  would  produce  a  specified  voltage  change  on 
the  monitor.   At  this  point  the  system  was  ready  to  run  and 
the  next  step  was  to  select  the  linearized  operational  mode 
on  the  monitor,  position  the  probe  in  the  free  stream,  and 
maximize  the  frequency  response  of  the  wire  by  setting  the 
inductance  to  the  minimum  value  that  would  not  cause  the 
signal  to  oscillate.   The  stability  of  the  system  was 
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monitored  during  this  step  on  a  dual-beam  oscilloscope  to 
ensure  no  oscillations  were  present  in  the  signal.   When 
the  probe  had  been  adjusted  to  the  flow  a  digital  voltmeter 
was  used  to  check  the  multiplier-averager  for  correct  opera- 
tion.  Then  the  x-y  plotter  was  set  to  record  the  signal  on 
an  appropriate  scale. 

Now  that  the  system  had  been  checked  for  correct  operation, 
it  was  positioned  as  close  to  the  upper  wall  as  practicable 
and  the  pen  was  dropped  on  the  plotter.   The  sensor  was  then 
driven  across  the  wake  to  a  point  near  the  bottom  wall,  and 
the  pen  lifted.   After  changing  the  paper  in  the  plotter, 
the  same  procedure  was  used  to  make  a  traverse  back  to  the 
top  of  the  test  section.   A  total  of  four  traverses  were 
made  at  each  station  to  record  both  the  mean  flow  and  the 
mean  square  of  the  fluctuating  flow.   Each  of  these  quan- 
tities was  recorded  at  two  different  scales  in  order  to 
facilitate  translation  into  digital  data  for  the  IBM-360. 

D.   TWO-DIMENSIONAL  ANEMOMETER  DATA 

The  basic  operation  and  accuracy  checks  of  the  x-array 
hot  wire  are  very  similar  to  those  of  the  single  normal 
wire.   However,  some  additional  complications  are  created 
by  having  two  separate  wires  on  the  probe.   The  terms  are 
now  defined, 

e  =  e  +  e' 
u  =  u  +  u' 
v  =  v  +  v' 
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As  explained  in  Appendix  B,  the  outputs  of  the  two  wires 
(A  &  B)  were 


eA  -  CA  (u  +  v) 


eB  =  CB  (u  -  v) 


Thus,  for  a  given  flow,  the  values  of  C.  and  Cn  cannot  be 

A  D 

determined  easily.   In  order  to  find  these  constants  one 
must  have 


CA  "  CB  "  C 


In  this  study  the  author  used  a  basic  physical  character- 
istic of  pipe  flow  to  set  these  constants  equal,  i.e.,  that 
the  lateral  velocity  must  be  zero  at  the  centerline.   Recall 
that,  at  y  =  0 


eA  =  V  =  CBU  =  eB 


Thus  the  calue  of  (e.-eD)  was  set  equal  to  zero  so  the 

A    D 

constant  (C)  could  be  defined  as  in  (III-5).   Note  that  this 
procedure  is  equivalent  to  setting  V  =  0  at  the  centerline 
since. 


efl  -  eD  =  C(u  +  v)  -  C(u  -  v)  =  2Cv 

A       D 


2H 


eA  +  eR  =  C(u  +  v)  +  c(u  -  v)  =  2Cu 


The  actual  position  of  the  centerline  was  determined  from 
the  pitot-static  data  taken  in  phase  1. 

Again,  the  wires  were  adjusted  to  the  individual  flov;  at 
each  station  and  the  data  were  taken  in  the  same  manner  as 
described  in  the  previous  section.   However,  several  more 
traverses  were  required  to  get  all  of  the  data  available 
to  the  x-array  hot  wire.   The  settings  on  the  correlator, 
and  the  corresponding  flow  parameters,  are  listed  below. 
All  barred  quantities  were  taken  with  the  correlator  on 
DC,  and  averaged  before  plotting.   The  primed  quantities 
were  taken  with  the  correlator  on  AC,  and  squared,  averaged, 
and  then  plotted. 


SETTING  RESULT 


(A+B)  u 


! 


(A-B)  v 


u+v 


B  u-v 


SETTING 

RESULT 

(A+B)'2 

u-2 

p 

(A-B)'2 

v'd 

A'2 

(u+v)'2 

B'2 

(u-v)'2 
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IV.   DISCUSSION  OF  RESULTS 

Since  this  study  was  intended  to  carry  out  an  investi- 
gation of  free  turbulent  mixing  in  a  constant  pressure 
gradient,  it  would  be  appropriate  to  discuss  how  close  the 
actual  experimental  environment  approached  that  desired  for 
the  study.   It  can  be  easily  seen  in  Fig.  21  that  a  wake 
velocity  defect  was  initially  created  by  the  cylinder  as 
desired,  and  the  magnitude  of  u'   in  Fig.  4l  clearly  shows 
that  the  flow  was,  in  fact,  turbulent.   However,  it  can 
also  be  seen  that  the  adverse  pressure  gradient  encourages 
the  formation  of  a  large  boundary  layer  which,  as  it  mixes 
downstream,  interacts  with  the  wake  defect,  tending  to 
raise  some  doubt  that  free  turbulent  mixing  is  still  present 
downstream  of  Station  6.   It  is  especially  significant  to 
note  that  the  turbulence  near  the  wall  builds  up  from  a 
negligible  value  at  Station  1  (see  Fig.  41-50)  to  a  value 
greater  than  that  induced  into  the  wake  by  the  cylinder 
midway  down  the  test  section.   As  expected,  the  wake  defect 
eventually  begins  to  mix  out  in  the  favorable  pressure 
gradient. 

By  comparing  the  stress  calculated  from  the  pitot-static 
data  with  that  from  the  hot-wire  anemometer,  one  can  easily 
see  that  in  both  cases  the  stress  is  zero  on  the  centerline, 
then  builds  up  to  a  maximum  as  the  velocity  gradient  nears 
the  region  of  its  maximum.   At  the  later  stations,  both 
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types  of  data  show  a  shear  stress  sign  reversal,  and  its 
magnitude  begins  to  increase  as  the  traverse  approaches  the 
edge  of  the  test  section.   This  trend  in  the  data  was  prob- 
ably reflecting  the  increasing  interaction  with  the  boundary 
layer  as  the  air  flowed  downstream. 

Due  to  the  greater  capabilities  of  the  hot-wire  probe, 
it  was  used  to  take  data  very  near  the  wall,  and  the  velocity 
"can  be  seen  to  fall  rapidly  in  the  vicinity  of  the  boundary 
layer.   The  hot-wire  data  (Figures  1-80)  were  then  used  to 
prepare  Table  I  to  display  the  value  of  the  correlation 
coefficient  (C   )  and  n  in  the  turbulent  flow  field.   Note 

that  the  correlation  coefficient  C   approaches  zero  at  the 

uv  1  ^ 

centerline,  and  reaches  its  maximum  where  the  shear  is  a 
maximum  as  predicted  by  Schlicting  (See  Fig.  18.^  Ref.  2). 

Considering  the  relation  between  kinetic  energy  and  the 
shear  stress,  it  is  easy  to  see  in  Table  II  that  both  param- 
eters reach  a  maximum  at  nearly  the  same  point  in  the  flow. 
Comparing  the  anemometer  with  the  pitot-static  data  in 
Appendix  F,  the  maximum  shear  there  also  occurs  at  the 
lateral  position  corresponding  to  the  kinetic  energy  maxi- 
mum.  Thus,  it  is  clear  that  at  least  some  significant  func- 
tional relationship  exists  between  them.   The  ideal  relation- 
ship would  be  that  postulated  by  the  theory  of  eddy  viscosity, 
that  the  two  parameters  differ  only  by  a  constant  that  is 
a  function  of  x.   The  slight  displacement  of  the  two  maxi- 
mums in  this  case  would  imply  that  the  relationship  between 
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k  and  C   is,  in  reality,  a  more  complex  function.   As  an 
additional  attempt  to  clarify  the  relation  between  kinetic 
energy  and  shear  stress,  equation  (7)  was  approximated  by 
taking  the  characteristic  velocity,  as  the  square  root  of 
the  kinetic  energy,  and 


S  =  C1  I   =  S(x)  (IV-1) 


As  a  result,  the  relation  for  Reynolds  stress  becomes, 


Ti       ^u    of    \    i  1/2  8u  ,_„  „v 

-k  =  yt  87  =  s(x)  k   87        <IV-2> 


The  values  for  equation  (IV- 2)  were  calculated  and  plotted 
versus  traverse  position  for  several  stations  from  the  data 
in  Table  II, which  was  derived  from  Figures  21-80.   The 
results  of  these  plots  are  shown  in  Figures  93-97.   Once 
again,  there  is  a  definite  similarity  between  the  results 
of  the  approximating  equation,  (7),  and  the  calculation  of 


u'v1  from  the  hot  wire  data.   In  this  case  the  shapes  of 
the  curves  were  exceptionally  similar,  but  the  maximums 
were  still  displaced  by  a  small  amount,  preventing  direct 
correlation  by  multiplication  by  a  constant.   The  author 
feels  that  a  more  comprehensive  analysis  of  the  data  in 
this  report  is  warranted  by  the  obvious  visual  relationship 
of  the  curves.   An  in-depth  analysis  might  well  uncover 
a  significantly  better  model  for  calculating  turbulent 
flow  fields  than  are  currently  in  use. 
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TABLE  I 
Calculation  of  The  Correlation 
Coefficient  and  lita 


Oa, 


e:-e: 


**   fe+o^i^-Q 


afer-r-e^ ) 


(eA+e&)'2  4-  {jeZ-ej* 


$tation  1 

Pnsi  +",1  nn 

C)                      i 

1 

10.5 

1 

-.6529 

\ 
1.310 

11.0 

-.1-881 

1.3lll 

ii.5 

-.14-353 

l.k9k 

12.0 

-.2092 

1.714-5 

12.5 

-o0319 

1.812 

13.0 

.1501 

1.782 

13.5 

.2857 

1.1*89 

H|..0 

.3112 

1.292 

iii.5 

.1631 

I.2I4.3 

Station   2 

9.5 

-.W5 

.9032 

10.0 

-.i|170 

1.250 

10.5 

-.2029 

1.012 

11.0 

-.2353 

.9706 

n.5 

-.2^70 

1.058 

12.0 

-.3108 

1.21+3 
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Position 


/ 


12.5 
13.0 

13.5 
II4.O 

111..  5 

15.0 
15.5 


-.0381 
.1172 

.1880 
.2576 
.0851 
.0683 

-.2887 


v. 
1.219 

1.2*49 

1.188 

1.31+3 
1.277 
1.153 
.875- 


9.0 
10.0 
11.0 
12.0 

13.0 
llf.O 

15.0 


Station   3 


-.31+01). 
-.3661 
-.2I4OI 

.0915 
.1889 
.1097 


.878 
.922 
1.122 
1.227 
1.281 
1.211 
1.071 


Station  k 


8.0 
9.0 
10.0 
11.0 
12.0 
13.0 

II4.O 


-.6876 

-.5699 
-.1*562 

-.3599 

-.1709 

.0983 

.2615 


1.108 
1.281 
1,361 
1.24-26 
1.2472 
1.14714 
1.14-85 


30 


ll_    (cont.) 



Position 

C 

n 

V 

15.0 

.2815 

1.31U 

16.0 

.214.00 

1.1914- 

17.0 

-.OI4.23 

1.036 

Station  5 

8.0 

-.2750 

I.O83 

9.0 

-.331+0 

1.176 

10.0 

-.3U70 

1.317 

11.0 

-.301+0 

1.3la 

12.0 

-.1508 

1.388 

13.0 

.0582 

U|27 

II4..O 

.1907 

1.390 

15.0 

.258ii. 

1.317 

16.0 

.18!].8 

1.161 

17.0 

-.0I4.18 

1.105 

Station  6 

7.0 

-.01514. 

i.ii+5 

8.0 

-.2951 

1.125 

9.0 

-.2809 

1.165 

10.0 

-.2711 

1.269 

11.0 

-.14.096 

1.119 

12.0 

-.0900 

1.2514. 

13.0 

.0295 

1.3714. 
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6  ( c  on  t ) 


Position 


ll|..0 
15.0 
16.0 
17.0 


.2051 
.2327 
.1701; 
.0)4.61 


Station  7 


1.296 
1.283 
1.212 
1.0l|.9 


7.0 
8.0 
9.0 
10.0 
11.0 
12.0 
13.0 
11+.0 
15.0 
16.0 
17.0 


.0681]. 

-.131*1 

-.2719 

-.3906 
-.2718 

-.15  A 

.0857 
.2013 
.1814.5 

.1131 
.6214.6 


1.103 
1.109 
1.158 
1.291 
1.279 
1.228 
1.2614. 
1.306 

I.2I4.3 

1.176 
I.OI4.8 
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Station  8 

Position 

C 

A 

7.0 

.1517 

\ 
l.ll+U 

8.0 

-.1027 

1.010 

9.0 

-.3392 

1.138 

10.0 

-.7895 

1.095 

11.0 

-.7298 

1.191*. 

12.0 

-.3909 

1.205 

13.0 

.1279 

1.165 

Ilj-.O 

•U5io 

1.232 

15.0 

.26514- 

1.200 

16.0 

.0631+ 

1.151 

17.0 

-.I4.36O 

1.185 

Station  10 

9.0 

.okSk 

1.095 

10.0 

-.7389 

1.130 

11.0 

-.81+39     • 

1.079 

12.0 

-.3313 

1.080 

13.0 

.2910 

1.091 

1U..0 

.l|8l8 

1.105 

15.0 

-.0229 

1.1U5 

16.0 

1  - 

-.0116 

1.020 
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Station  9 


Position 


t 


9.0 
10.0 
11.0 
12.0 

13.0 

11*. 0 

15.0 


-.3611*- 
-.61*85 

-.7565 
-.3710 

.1835 
.5671* 
.3130 


I.123 
I.063 
1.130 
1.091 

1.11*3 
1.11+6 
1.11*0 
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TABLE   II 
Calculation  of  Kinetic  Energy 
And  Shear  Stress  Parameters 


uV  = 


Pfei; 


e 


:y 


'■-M^f^X] 


Station   1 

Position 

k 

u'v» 

12.8 

993.9 

36,81 

13.0 

96^.0 

70.99 

13.2 

928.2 

8[|..1)| 

13.U 

870.3 

9U.66 

13.6 

780.9 

115.69 

13. 8 

657.3 

89.^0 

II4..0 

512.7 

78. 8S 

1I-.2 

331.3 

U7.33 

Ik.k 

239.3 

A.  18 

34.6 

155.1 

.  13.15 

Station  2 

12.0 

535.0 

-83.O 

12.5 

5*4-5.0 


10.i| 
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Station   2    (cont 

..) 

Position 

k 

u'v' 

12.8 

534-0 

20.74 

13.0 

531.-'+ 

31.11 

13.2 

521.0 

36.29 

13.4 

505.5 

31.11 

13.6 

L166.6 

20.74 

13.8 

414.8 

7.78 

1ft..  0 

362.9 

5.18 

1J+.2 

3H.1 

2.59 

14  -k 

267.0 

5.18 

14.6 

223.0 

5.18 

IJ4..8 

189.3 

7.78 

15.0 

153.0 

5.18 

Station  3 

12.  k 

378.5 

-22.08 

12.6 

376.9 

-11.04 

12.8 

376.9 

3.15 

13.0 

380.1 

15.77 

13.2 

381.7 

28.39 

13.4 

378.5 

28.39 

13.6 

364.3 

29.97 

13.8 

343.8 

29.97 
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Station   3    (cont.) 


Position 


k 


u'  v' 


ll+.O 
34.2 
14.  k 
14.6 

II4-.8 
15.0 
1^.2 
l5.Ii 


320.2 

296.5 
265.0 

236.6 
205.0 

176.6 
135.6 
119.9 


29.97 
23.66 
22.08 

17.35 

12.62 

9. 1+6 

6.31 

3.15 


Station  1+ 


13.0 
13.2 

13.U 

13.6 

13.8 
II]..  0 

11^.2 

li|.i+ 
34.6 
II4..8 
15.0 
15.2 


235.9 

237.5 
238.2 

237.5 
235.1 
232.8 

226.7 
219.0 
209.8 
199.0 
186.0 
172.9 


11.53 
16.  lit. 

21.52 
25.36 
27.66 

29.97 
32.28 
30. 7k 
25.36 

23.05 
19.98 

16.91 
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■■  ■                              — 

Station  4    (cont 

>.) 

Position 

k 

u'v' 

I5.li 

156.8 

16.91 

15.6 

139.9 

12.30 

15.8 

no    i 

10.76 

16.0 

101.  4 

6. 45 

16.2 

86.1 

7.69 

16.1] 

72.2 

6.15 

16.6 

61.5 

3.84 

16.8 

52.3 

3.84 

Station  5 

12.0 

217.5 

-14.32 

12.2 

214.8 

-   8.05 

12.1+ 

213.0 

-  3.58 

12.6 

213.0 

0.90 

12.8 

215.7 

4.48 

13.0 

217.5 

7.16 

13.2 

221.0 

II.63 

13.1+ 

224.6 

14.32 

13.6 

228.2 

17.00 

13.8 

230.0 

19.69 

1U.0 

232.7 

21.43 

14.2 

231.8 

24.16 

li+.lt- 

,     229.1 

25.95 
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Station   5    (cont 

.) 

Position 

k 

u'v' 

III.  6 

226.1* 

— "  -  —  —■—-■—  —  ——■   ■■■  — ..-.—.    .,._.__-. 

25.95 

UwB 

221.9 

25.95 

15.0 

21)].  8 

26.85 

15.2 

208.5 

26.85 

15.4 

197.8 

26.85 

15.6 

187.9 

2i|.l6 

15.8 

176.3 

21.1|8 

16.0 

166.5 

17.90 

16.2 

152.1 

14.32 

16.  k 

137.8 

H.63 

16.6 

125.3 

8.95 

16.8 

107.  k 

1.79 

17.0 

93.97 

-   1.79 

17.2 

85.0 

-   2.69 

17.14- 

73.8 

-   3.58 

Station   6 

•- 

13.0 

191.1 

2.71 

13.2 

197.5 

6.3)4- 

13.il 

202.0 

10.87 

13.6 

205.6 

13.59 

13.8 

203.  k 

18.12 

1U.0 

211.1 

20.81+ 
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Station  6  ( coat. ) 


Position 


k 


1  vt 


U'  V 


11+.2 

lli-.ll. 
H4..6 
14.8 

15.0 
15.2 

l$.k 

15.6 

15.8 

16.0 
16.2 
16.  k 
16.6 
16.6 
17.0 
17.2 
17.4 


212.0 
212.0 
211.1 
208 .1; 
201+.7 

199.3 
193.9 
188.1* 
181.2 
171+.8 
16b  .0 
153.1 

135.0 

129.5 
121.1). 

116.9 


Station  7 


23.55 

2]|.!46 

23.55 
22.65 

21.71] 

22.65 
19.02 
16.31 
13.59 
10.87 
7.25 

5.144 
3.62 
2.72 
0.0 
•   3.62 


13.0 
13.2 

13.1* 

13.6 


169.3 
173.9 
178.1+ 

182.9 


10.37 

111.  ^9 
17.21 
18.11 
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Station   7    (cont 

.) 

Position 

k 

u'v' 

13.8 

188.  b, 

18.11 

11*.  0 

192.9 

18.11 

111.  2 

195.6 

18.11 

14.4 

193.3 

18.11 

14.6 

198.3 

18.11 

34.8 

215.0 

18.11 

i5.o 

193.8 

17.21 

15.2 

190.2 

17.21 

15.4 

184.7 

15.1)0 

15.6 

180.2 

Hu49 

15.8 

173.9 

11.77 

16.0 

169.3 

9.06 

16.2 

163.9 

8.15 

16.1* 

158.5 

7.25 

16.6 

154.0 

6.34 

16.8 

151.3 

4.53 

17.0 

148.5 

1.81 

17.2 

11+7.6 

-   3.62 

17.1+ 

11+9.1+ 

-10.87 

17.6 

154.9 

-16.30 
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Station  8 

Position 

k 

U  '  V  ' 

15.0 

171.8 

11.31* 

15.2 

170.1 

9.60 

15.1* 

166.6 

9.60 

15.6 

161*.  0 

7.85 

15.8 

159.6 

5.23 

16.0 

156.1 

2.62 

16.2 

151.8 

-     1.75 

16.  k 

11+8.3 

-     5.23 

16.6 

ll4-7.li- 

-      9.60 

16.6 

150.0 

-   111. 83 

17.0 

160.5 

-  18.32 

Station  9 

12.0 

136.  k 

-13.22 

12.2 

133.1 

-   9.92 

12.  k 

131.1* 

-  5.79 

12.6 

131-U 

-   1.65 

12.8 

133.9 

2.1*8 

13.0 

138.8 

6.61 

13.2 

114-0.5 

10.71* 

13.U 

11*3.8 

11*. 05 

13.6 

11*7.1 

17.36 

13.8 

151.2 

20.66 
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Station   9    (cont 

.) 

Position 

k 

u'v' 

ll+.O 

152.9 

21.49 

34.2 

152.1 

20.66 

14-4 

ll(9.6 

19.83 

II4..6 

147.1 

16.53 

14.8 

143.8 

14.88 

15.0 

ll|2.1 

11.57 

Station   10 

12.0 

123.0 

-10.66 

12.2 

119.7 

-   6.56 

12.  k 

119.7 

-  1.64 

12.6 

121.4 

2.46 

12.8 

124.6 

6.56 

13.0 

126.3 

9.84 

13.2 

132.0 

13.914 

13.  k 

135.3 

16.40 

13.6 

138.6 

18.04 

13.8 

l4i.i 

19.68 

II4..O 

140.2 

17.22 

14.2 

138.6 

15.58 

ll4-.it- 

138.6 

12.30 

14.6 

137.8 

7.38 

14.8 

136.1 

3.28 
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FIGURE   Si 

The  wind  Tunnel  Inlet 


12  4 


FIGURE  82 


The  Wind  Tunnel  Test  Section 
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P]  3UHE   83 


The  wind  Tunnel  31cwer 
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PIGURi  QU. 


The  Test  Section  Traverse 
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FIGURE   85 

The  Test  Section  Access  Panels  and 
Static  Pressure  Ports 
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The  Manometer  3oard" 
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FIGURE   87 

The  Wake  Generator 
With  Abrasive  Strip 
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FIGURE   88 


The  Data  Sensors 
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F I  (JURE   89 


The  Test  Section  During  Phase  I 


132 


FIGURE   90 

X-Array  and  formal  Wire 
Anemometers 
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FIGURE   91 

The  Anemometer  Modules 
and  Correlator 
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The  Test  Section  Daring 


Phases  II  and  III 
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FIGURE   93 
Station  2  Turbulent  Stress  Comoarison 
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APPENDIX  A 
DEVELOPMENT  OF  THE  GOVERNING  EQUATIONS 
OF  MOTION  FOR  TURBULENT  INCOMPRESSIBLE  FLOW 

First,  the  following  equations  are  the  basic  Navier- 
Stokes  equation  for  incompressible  Newtonian  flow: 
continuity : 

V  •  V  =  0  (A-l) 

momentum: 


■°L   +  v  •  V  V  =  -  -  VP  +  v  V      V        (A-2) 
ou  p 


Expanding  the  DEL  terms  they  are 


|u  +  |S  +  |*  =  o  (A-3) 

3x    3y    3z 


2      2      2 
3u  ,   3u    3u  ,   3u    1  3P  ,   ,8  u  .  3  u  ,  9  U\  /«  i,  \ 
sr  +  u-s—  +  v^—  +  w-s— ■  =  --  ^—  +  v( — 5-  +  — K-  +  — x-)  (A-4) 
at     3x     3y     3z     p  3x      9x2    9y2    ^2 


3v     3v     3v     3v     1  3P  ,   ,32v    32v  ,  32Vv  ,.  0 
3t     3x    3y     3z     p  3y     ^2        9y2    8z2 


2      2  2 

3w    3w    3w    3w  =  —  -jr—  +  v( — ~-  +  — ~  +  — 5-;  CA-o; 

"5*r  +  Utt—  +  Vtt—  +  Wt-    p  3z      ~  2    a  2  ~    2 

3t     3x     3y     3z  3x     3y  3z 
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Now  substitute  equation  (A-l)  into  (A-1!),  (A-5),  and  (A-6) 
to  get 


3u  9   2  ,  3  „„  ,  3  „TT  _   1  3P  ,   ,32u  ,  32u  ,  32wN 

3t  +  3xU   +  3yUV  +  *zuw  "  "p  37  +  ^(~2  +  ~2  +  ~2 > 

dX  dy  dZ 

__J  +  3  uv  +  i  v2    3  _,  J.  3P     (3fv  +  9^u  1fv) 

at  dx      ByV     3z™     p  3y    *<ax2  8y2  ^ 

8w  ,  3       8  _„_    3   2     1  3P     ,32w  ,  32w  ,  32w^ 

_  +  _xUW  +  _yVW  +  _zW  _____  +  v(_  +       +  __) 


With  p  and  v  constant  and  using  the  following  definitions 

u=u+u'  w=w+w' 

v=v+v'  P=P+PT 

equation  (A-3)  becomes 


|v(u  +  u')  +  f-  (v  +  v')  +  |-_(w  +  w«)  -  0     (A-7) 

dX  dy  d  Z 


Since  the  integral  over  time  of  any  fluctuating  quantity  is 
zero,  taking  the  time  average  of  equation  (A-7)  results  in 


3u    3v    3w  (A_8) 

3x    3y    3z 


Then  subtract  this  from  (A-7)  to  get 


j_!  +  |_1  +  |s_  _  o  (a-9) 

3x    3y     3z 
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Now  substitute  the  definitions  above  into  equations  (A-4), 
then 


f-.Cu+u')  +  |-  (u+u')2  +  f-  (u+u»)(v+v')  +  |-(u+u!)(w+w') 

do  dX  a j  dZ 

2  2  2 

=  4i(P+p,)  +  vC— 2(u+u')  +  i_2(u+u')  +  i-  (u+u«)] 
p  9x  9y  9z 

(A-10) 
and  after  taking  the  time  average 


!£+-|u"2  +  !-u'2+!-uT+jL  HTTf"  -r  |  ot  +  |  u'sT1"  (A-ll) 
8t    9x      9x       9y      9y        9z      9z 


_   1  9P  ,   r92u  .  92u  ,  32u-, 
H         9x     9y     az 


In  a  similar  manner,  one  can  get 


I?  +  !■  uv  +  !■  u^v7"  +  |  v2  +  I  v'2  +  i  vw  +  I-   v'w'  (A-12) 
9t    9x     9x       9y     9y      9z     9z 

2—    2—    2— 
_   1  9P  ,   r9  v    9  v  ,  9  vn 

-  — -  -r—    +  VL q  +  q  +  oJ 

p  8y      9x2    9y2    9z2 


If  +  tx™  +  fx^  +  fy™  +  lyv'w'  +  !zw2  +  |zw'2  (A~13) 

_   1  9P     r92w    92w    92w-, 
p  3Z      9X2    9y2    9z2 


In  the  special  case  for  a  two-dimensional  mean  flow  with  a 
steady  mean,  equations  (A-ll)  and  (A-12)  can  be  written 


143 


3u2    A    8    — -  „      18    /-x...i2x         9    ,rr7T7TN     ,      .r92u    ,    92i 


+      7T     UV 


1    9    /-~.       1 2 >.         9    / — j — j-v     .       r9    u    .    9    u-i 


9x  9y"v  p    3x^'^      '         9yv"    v    '     '     VL„2    '    .2 

ox  dy 


(A-14) 


9uv    ,     9      2  13    ,-        ,2v         9    (— = — r>  r92v        92vn 

97"  +  Tyv     =  "p  9"y(p+pv     }   "  ¥x(u  v,)  +  v[7-2  +  7T] 

d  X     d  y 

(A-15) 

Then  applying  the  standard  boundary  layer  assumptions 

P  =  P(x)      from   |=-PUe^  (A-16) 

and  in  addition  make  the  following  approximations 

p  u,2<<  P  (A-17) 


2-    2- 
v[i_U  +  3U-J  ^  1  ulyT  (A-18) 

9x^   9y^     9y 


When  these  assumptions  are  substituted  into  (A-l'O,  the 
following  equation  results: 


2  "" 

oil    .  a    J.  d  r     9  — i f-  /  «  -,  p.  \ 

33T  +  3yuv  =  "p  3?  "  Tyu  v  (A-19) 


Which  can  be  manipulated  into  the  form  in  equation  (A-20) 
using  continuity 


„  |a  +  7  |a  =  -i  |£  -  »  u^  (A-20 

9x     9y    p   9x   9y 
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applying  the  Boussinesq  laminarization,  the  turbulent  shear 
becomes 


fy"'V  =^Ty  =  I(pt(x)  ff)  CA-21) 


then,  x  =  -p  iT'T7",  is  the  Reynolds   Stress,  and  y,(x),  the 
Eddy  Viscosity,  is  a  function  of  the  flow  field.   If  (A-16) 
and  (A-21)  are  substituted  into  (A-20)  the  result  is 


srr     sT      dU  .2- 

"~  °u  .  —  dU    TT     e  ,     /  s  3  U  /»  ~~\ 

u  87  +  v  87  =  Ue  OT-  +  VT(x)  77  (A"22) 

J  8y 


An  alternate  approach  to  solving  these  equations  is 
possible  without  the  boundary  layer  assumptions  (A-l6,  17, 
&  18).   Since  there  are  three  equations,  (A-8),  (A-lM) ,  and 


_  _  _    2     2      

(A-15)  and  six  unknowns  u,  v,  p,  u'  ,  v'  ,  and  u'v',  three 

additional  equations  are  required.   By  subtracting  the  mean 

equations  from  the  full  equations  as  shown  in  (A-10),  three 

additional  momentum  equations  for  the  fluctuating  portion 

of  the  flow  result;  however,  several  additional  unknowns 

also  appear: 


"5-5  2 

u'  ,   v'  ,  and  u'v'  . 


One  might  also  formulate  additional  equations  for  these 
variables,  but  higher  order  correlation  unknowns  would 
appear  in  them,  and  additional  equations  would  be  required 
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ad  infinitum.   It  is  customary  at  this  point  then  to  approxi- 


,2    ,2 


mate  v'  ,  u'  ,  and  u'v'  by  some  function  of  the  mean  flow 
and  solve  the  system  of  equations. 
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APPENDIX  B 

DEVELOPMENT  OF  THE  EQUATIONS  FOR 
THE  HOT-WIRE  ANEMOMETER 

I.   Single  Normal  Hot-Wire 

First  assume  the  following  definitions 
e  =  e  +  e ' 
u  =  u  +  u1 

Then,  since  the  linearized  output  of  the  hot-wire  is 
proportional  to  u 

e  =  C  u  (B-l) 

or 

e  +  ef  =  C(u  +  u' )  (B-2) 

Setting  averages  equal  to  averages,  since  they  must  be 
independent  of  the  fluctuations, 

e  =  C  u  (B-3) 

or  solving  for  u 


u  =  1  e  =  ^  e  (B-H) 

ee 
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and 


U_  =  e_ 

U     e 


Now  setting  the  fluctuations  equal 


>  =  r    11 1 


u'  =  i  e' 
C 


,2    1    ,2 

u '   =  —  ?  e  ' 

C 


which  from  (B-4)  becomes 


U  2 


e 


e 


or 


(B-5) 


C  u'  (B-6) 


uft  =  -^2  e'c  (B-7) 


1 2     ,2 

u'   _  e1 

—2    — 2 
u     e 
e     e 


II.   X-array  Hot-Wire 

Now  assume  in  addition,  that 


v  =  v  +  V      eA  =  eA  +  eA»    eB  =  eB  +  e- 


(B-8) 


(B-9) 
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then  since  the  wires  are  separated  by  90° 


¥  usin0  +  vcos0 


ucosG  +  vsin0 


wire  A 


wire  B 


Now  note  that  6  =  ^°  so 


eA  =  CA  (  U  +  v) 


eB  =  CB  (  u  -  v) 


(B-10) 


At  the  centerline  of  the  flow  the  lateral  velocity  is 
approximately  zero,  so  e.  =  eR  and,  therefore,  C.  =  CR  -  C. 
For  this  reason,  once  the  probe  is  adjusted  to  force  e.  =  eB 
at  the  centerline  the  following  relations  hold 


eA  =  C  (u  +  v) 


eB  =  C  (u  -  v) 


(B-ll) 
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But  if  these  equations  are  summed  and  solved  for  u 


u  =  ~  [eA  +  eB]  (B-12) 


and  similarly 


v  =  k   [eA  "  eB]  (B"13) 


Then  if  3  =  1/2  C,  and  the  definitions  of  the  parameters  are 
substituted  into  (B-12)  &  (B-13),  the  resulting  equations 
are 


u  +  u'  =  3[eA  +  eA'  +  eg  +  eB' ]  (B-14) 


v  +  v'  =  B[eA  +  e  '  -  eR  -  eBf]  (B-15) 


As  in  B-I,  set  averages  equal  to  averages  and 


u  =  3[eA  +  eB]  =  3[eA  +  eB]  (B-16) 


v  =  3[eA  -  eB]  =  B[eA  -  efi]  (B-17) 


setting  fluctuations  equal  to  fluctuations, 


u'  =  BCeA'  +  eB']  =  3[(eA  +  eB)f]  (B-l8) 


v'  =  3[eA»  -  eB']  =  3[(eA  -  eB)*]  (B-19) 
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Then  square  and  average  to  get 


u'2  =  32[(eA  +  eB)'2]  (B-20) 


v*2  =  32[(eA  -  eB)'2]  (B-21) 


Returning  to  equations  (B-ll)  and  setting  primes  equal  to 
primes 


eA'  =  C(u'  +  v' ) 


eB«  =  c(u»  -  V) 


square  and  average 


eA'2  =  C2  (u'2  +  2u»v'  +  v'2) 


e  '2  =  C2  (u'2  -  2u'v'  +  v'2) 

D 


then  subtract 


or 


Now  add  equations  (B-23) 
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(B-22) 


(B-23) 


e'2  -  e  '2  =  4C2u'v'  (B-24) 

A  D 


^  [eA-2  -  eB'2]  (B-25) 


and 


eA'2  +  eB'2  =  2C2[u'2  +  v'2] 


or 


u'2  +  v'2  =  232[eA'2  +  eB'2]  (B-23) 


Then  sum  (B-20)  and  (B-21) 


u«2  +  v«2  =  32[(eA  +  eB)'2  +  (eA  -  eB)'2]    (B-24) 


equating  (B-23)  and  (B-24),  and  rearranging 


2[e  «2  +  e R'2] 
A 5 =  1  (B-25) 


(eA  +  eB)'2  +  (eA  -  eB)<2 


Where 


PTe  "  Pe  "  l/2P["e2  +  ve^ 


or 


Pm  -  P 


2[-lS— -£]  =  (l^)2  (e  ♦  eB)2  +  (-L)2  (e  -  eB>2 


SO 


p    -  p 
32  =    2(rTe    e) 


p[(eA  +  eB)2  +■  (eA  -  efi)2]e 


T  r-. 


and 


2(PTe  "  Pe) 


is  a  known  array  of  data  from  the  pitot-static  phase  of  the 
study. 
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APPENDIX  C 
DEVELOPMENT  OF  AN  INTEGRATED  FORM  OF  THE 
MOMENTUM  EQUATION  FOR  THE  CALCULATION  OF 
THE  LOCAL  TURBULENT  SHEAR  COEFFICIENT 

First,  the  basic  equations  of  motion  for  a  two-dimensional 
mean  flow  are 

continuity : 


u   +  v   =  0  (C-l) 

x    y 


momentum: 


uu   +  vu   =  -  -  P   +  -   t  (C-2) 

x     y     p  x   p  y 


then  integrate  equation  (C-l)  from  zero  to  a  point  y 


u  =  -  v 
x      y 


/  u   dy  =  -/   v  dy 
0  0   y 


y  y 

/  uY  dy  =  -  [v  +  f(x)] 
0   x  0 


thus 


y 

v(y)  =  -  /  u  dy  +  v 

0  «y" 


15*J 


But  v     is  zero,  so 


y 
v(y)  =  -  /  "  uY  dy  (C-3) 

0   x 


Now  simply  and  multiply  equation  (C-l)  by  u  and  add  to 
(C-2)  to  get 

2uu   +  (vu   +  uv  )  =  -  -  P   +  -^  (C~il) 

x     y    y     p  x   P 

where  the  terms  on  the  left  can  be  recognized  as  differentials 
of  the  form 


2 
2uux  =  (u  )x 


vu   +  uv   =  (uv) 

y    y     y 


thus  equation  (C-*0  becomes 


,      x 

2uu  +  (uv)   =  -  -  P  +  -2- 
x      yy     p   x    p 

and,  since  v(0)  =  0,  when  this  equation  is  integrated  from 
y  =  0  to  an  arbitrary  point  y  =  y  , 

y         y         y 

/   2uuv  dy  +  /  T(uv)v  dy  =  /  T2uu   dy  +  u(y  )v(y  ) 
0  0      y      0 


and  the  result,  when  equation  (C-3)  is  substituted  for 
v(y  ) ,  is  the  following  relation 
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yT  YT  YtP  T(y     ) 

/      2uu   dy    -    u(y    )    /      u   dy   -   -   /      -^dy    +   — -I- 
0  T0  Op  p 


Which   can  be   rearranged   into   the    following   form, 

yT  yT 

/      u    [2u-u(y    )]dy   -   -  I     /      P   dy   +  I  T(y    )    (c-5) 
0X  T  pQxpT 


And  can  be  solved  for  the  coefficient  of  shear  at  y 

T 
PU„ 


yT  yT 

CT(yT)  =  —V"  /   pxdy  +  [-V  /   u  [2u-u(y)]dy] 
T   T     pu/   0    x      u/   0     x       T 

e  e  (C-6) 


An  alternate  approach  to  this  development  would  be  to 
substitute  equation  (C-3)  into  (C-2) 


-  9u    du       /,3Uv,r  _    1  9P  .  1  9t  (r    t\ 

dx    9y  q    9x        p  8x    p  9y 


Then  non-dimensionalize  by  division  by  the  freestream 
velocity  squared,  and  use  the  product  rule  of  the  calculus 
to  get  the  following  form  of  the  momentum  equation: 


u  t   u  >     (   u  \2   Uex    (   u  v    >y ,   u.x  ,   /  u  v  /_exN  ryr  u  n  - 

u^iT^x    (u~}   u~ (~)    !    (u~}  dy-(~)  (u~")  ;  (~)dy 

e  ue  x     e   i  e      e  y  0    e  x     e  y  e   0   ue 

C-  (C-8) 


-^) 


2    x 
p  u  2 
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2t 
Where  C  =  p  >    anc*  tne  subscripts  x,  and  y  indicate 

T   P  u 
H   e 

partial  differentiation  with  respect  to  x,  and  y  respectively 
Now,  if  the  condition  that  equation  (C-8)  be  satisfied  at 
every  point  in  the  flow  is  relaxed  to  allow  it  to  be  satis- 
fied on  the  average  between  y  =  0  and  an  arbitrary  point 
y  ,  then  (C-8)  can  be  integrated  by  parts  to  get  the 
following  relation  for  C 


C  (y  )    yx       u(y  )  u 

_I — 1_  =  /   [2(-^-) —][(—)  +(iL)(_ex_)]dy 

vu    u      u     u   u 


2      0      ee      exee 

y^   1 

0   2pu^ 

Now  if  the  original  definition  of  the  shear  term 


(C-9) 


T  =  -  P  U'v' 


is  substituted  back  into  equations  (C-5)  and  (C-9),  the 
following  relations  result: 

unrv^"  =   /   u Y[u(y  )  -  2u]dy  -  ±  f      P   dy    (C-10) 
0    x    t  p  Q     x 

u-^  =  u  2  /  T[^l-  -  2(JL)][(JL)  +(JL)(^S£)]dy 

e   0     u         e     exee 
e 

yTP  cc-ii) 
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Finally,  to  check  the  validity  of  these  equations,  let 
y   go  to  infinity  in  the  limit  and  assume  that 


P   =  -  p  u  u 

e   ex 


Then,  equation  (C-5)  becomes 


2   /  uu  dy  -  u    /  u  dy  -   /  u   u  dy  =  0 
0    x      e  p    ?  0      ex 


Since   the   shear,    and  therefore   C    ,    is    zero   at   °°  and  then 

T 

CO  2  °°  °° 

2      /    (%-)      dy   -      /    [(u   u)      -   uup    ]dy   -      /   u   u      dy=0 
O^x  0  x  x  0eex 


or 


/    [(iT)      -    (u   u)      +   uu        -   u   u    ]   dy   =    0 
Ox  xx 


thus,  combining  terms 


CO  CO 

/  IX.   ~-  (~   1)]  dy  +   /  uquo  (Ji-   l)dy  =  0 

(C-12) 


e   u   u         "    n   e  e„  u 
0        e   e      x     0      x   e 


and  if  the  flowing  parameters  are  defined 


0  ue  ue  0   ue 
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then  equation  (C-12)  may  be  written 


[u  26]   +  u  u   6*  =  0 
e        e  e 
x       x 


which  when  the  first  term  is  differentiated  becomes 


u  2  J^-  +  2  u  u    6  +  uu   6*  =  0 
e   dx      e   e       e   e 

x  x 


2 
or  dividing  by  u 


u 
e. 


dx   u 
e 


[26  +  6  ]  =  0  (C-13) 


which  is  simply  the  well-known  momentum  integral  equation 
for  a  boundary  layer.   A  similar  procedure  will  reduce 
equation  (C-9)  to  (C-13)  as  well. 
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APPENDIX  D 

DATA  CHECKS 

In  an  experimental  study,  it  is  very  important  to  have 
some  idea  of  the  quality  of  the  data  that  results  from  the 
experiment.  In  this  Appendix,  several  methods  of  checking 
the  accuracy  of  the  data  in  this  study  are  discussed. 

CONTINUITY  CHECKS 

Using  the  continuity  equation,  it  is  apparent  that, 

since  the  mass  flow  must  be  constant, 

wall 
/    udy  =  constant  (D-l) 

0 

Therefore,  this  equation  was  added  to  the  program  described 
in  Appendix  E.   The  results  of  this  check  for  the  pitot- 
static  data  varied  as  much  as  30$  over  the  entire  test 
section,  yet,  if  the  largest  and  smallest  value  are  neglected 
the  range  is  narrowed  to  ±  10$.   The  author  believes  that 
much  of  this  difference  is  due  to  the  inability  to  position 
the  pitot-static  probe  near  the  wall  in  some  cases,  while 
in  others,  the  head  of  the  probe  could  be  placed  adjacent 
to  the  wall  surface.   This  problem  is  especially  critical 
at  Station  1  where  the  boundary  layer  is  thin,  and  mass 
flow  near  the  wall,  at  free  stream  velocity,  comprises  a 
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significant  part  of  the  total  mass  flow  in  the  tunnel 
See  the  computer  output  in  Appendix  P. 

MOMENTUM  CHECKS 

Recall  the  integral  momentum  equation  derived  in 
Appendix  C, 


u 
§£  +  — -  [29  +  6*]  =  0  (C-13) 

ue 


which  is  a  special  case  for  the.  free  stream  where  C  is 
zero.   Taking  this  equation  as 


El 

— 

de 

dx 

du 

E2 

= 

1 

e 

u 

dx 

e 

E3  =  2G  +  6 

then  the  non-dimensional  residual  of  this  equation  was 
programmed  as 

Epsilon  =  (El  +  E2  *  E3)/(E1  -  E2  *  E3)      (D-2) 

and  its  values  are  listed  in  the  computer  output. 

In  general,  this  ratio  showed  the  data  to  correlate 
fairly  well  with  the  restrictions  imposed  by  the  momentum 
integral  equation  with  an  average  residual  of  less  than 
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25%   and  as  little  as  +  15%   for  most  stations;  however, 
Station  7  data  was  completely  out  of  line  with  a  residual 
greater  than  two.   The  primary  reason  for  the  failure  of 
Station  7  to  check  satisfactorily  was  an  inclusion  of  a 
significant  portion  of  data  at  the  upper  portion  of  the 
traverse  that  reflected  boundary  layer  interference.   The 
inclusion  of  the  resulting  velocity  decay  in  the  calcula- 
tion of  the  momentum,  and  displacement  thicknesses  created 
a  large  error  in  their  values.   This  error  is  almost 
certainly  the  cause  of  slightly  high  values  of  Station  8 
as  well. 

ENERGY  CHECKS 

Having  calculated  the  local  turbulent  shear  coefficient, 
it  was  possible  to  also  integrate  to  the  residual  in  the 
energy  integral  equation, 


^(u  V)  =  2/1  (|H)dy  (D-3) 

dx  e  3  n   p   dy 


where  the  energy  thickness  <S  _  is  defined  as 


oo 

6,  ~      f   u(u  2  -  u^)dy  (D-4) 

^0     e 


Therefore  taking  the  term  on  the  left  as  A,  and  that  on 
the  right  as  B,  the  relation  was  programmed  and, 

ZETA  =  (A  -  B)/(A  +  B)  (D-4) 
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The  results  of  this  data  check  were  inconclusive.   The 
magnitude  of  one  of  two  terms  above  (A  and  B)  was  so  much 
larger  than  the  other  that  the  difference  over  the  sum  was 
due  solely  to  an  error  in  the  programming  vice  an  inaccuracy 
in  the  rav;  data. 

HOT-WIRE  CHECKS 

In  Appendix  B  it  was  shown  in  equation  (B-25)  that  these 
terms  were  set  equal  to  n>  and  values  for  n  are  listed  in 
Table  I.   It  was  immediately  obvious  in  this  case  that  the 
error,  although  not  large,  was  directly  related  to  the 
strength  of  the  turbulence  in  the  flow.   Thus,  the  magni- 
tude of  the  lack  of  algebraic  correlation  in  this  check 
may  have  been  due  to  an  inability  of  the  equipment  to 
handle  the  rapid  variations  caused  by  severe  turbulence. 
It  is  important  to  note  here,  that  a  change  of  as  little 
as  10$  in  the  individual  values  of  the  parameters  used 
to  calculate  n  would  bring  even  the  worst  value  of  n  back 
to  within  10JS  of  1.0.   The  author  felt  that  the  data  from 
Phase  III  of  the  experiment  was  well  within  the  bounds  of 
acceptable  experimental  data  spread. 
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APPENDIX  E 
DATA  REDUCTION  COMPUTER  PROGRAM 

A .  GENERAL 

The  Fortran  IV  computer  program  included  in  this  Appen- 
dix is  a  revised  form  of  that  developed  by  the  study  of 
Reference  1 .   The  program  inputs  the  pitot-static  data  of 
phase  I  of  the  experiment,  and  a  limited  amount  of  data 
from  phase  II  to  calculate  the  wake  momentum  thickness, 
displacement  thickness,  energy  thickness,  local  turbulent 
shear  coefficient,  and  Eddy  Viscosity. 

At  the  inception  of  the  present  study  the  portion  of 
the  program  dealing  with  the  turbulent  shear  coefficient 
did  not  produce  accurate  values  for  C  .   The  author  of 
Reference  1  was  unable  to  find  a  reason  for  this  error 
since  the  equations  appeared  correct,  and  there  were  no 
obvious  coding  errors.   He  postulated  that  the  unreasonable 
values  were  a  result  of  an  inability  of  a  forward  finite 
difference  scheme  to  accurately  approximate  the  partial 
derivatives  in  equation  (C-9)  of  Appendix  C.   Since  this 
appeared  to  be  a  reasonable  theory,  this  author  attempted 
to  obtain  values  for  the  derivatives  by  two  additional 
methods : 

1.   A  least  squares  polynomial  approximation  was  fitted 
to  the  data  at  a  given  position  y  from  the 
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centerline  then  this  polynomial  was  differentiated. 

2.   A  central  difference  scheme. 

Unfortunately,  these  improvements  resulted  in  no  signi- 
ficant change  in  the  values  of  t.   Thus,  the  values  for 
the  various  derivatives  were  printed  out  in  detail,  and 
an  order  of  magnitude  analysis  was  carried  out.   At  this 
point,  it  was  not  difficult  to  see  that  the  values  of  the 
parameter  U  were  exceptionally  smooth  in  the  x-direction, 
and  obviously  not  a  problem  to  approximate,  yet  the  values 
calculated  by  the  IBM-36O  were  not  even  similar  to  those 
calculated  by  hand.   The  author  then  attempted  to  follow 
the  instructions  programmed  for  the  computer  and  quickly 
discovered  that  a  subtle  error  was  hidden  in  the  program. 
Since  each  station  had  a  different  number  of  points,  the 
index  corresponding  to  a  point  y  above  the  centerline  was 
different  at  each  station;  therefore  the  arrays  requiring 
differentiation  were  read  into  a  dummy  variable  where 
J  =  1  corresponded  to  the  centerline  for  each  station. 
This  correction  produced  values  for  the  local  shear  coeffi- 
cient that  were  of  the  right  order  of  magnitude,  and  the 
central  difference  scheme  was  considered  to  produce  the 
best  values  for  the  derivative.   Some  of  the  hot-wire  data 
was  also  analyzed  using  this  program,  but  the  values  for 
C  were  not  very  similar  to  those  calculated  from  the  pitot- 
static  data.   Since  the  pressure  term  was  found  to  be  the 
driving  force  of  the  equation  for  C  during  the  order  of 
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magnitude  analysis,  the  author  felt  that  approximating 
the  static  pressure  data  with  the  data  from  phase  one  would 
be  unlikely  to  produce  realistic  values  for  the  parameters, 
and  no  additional  attempt  was  made  to  use  the  IBM-360  to 
reduce  the  hot-wire  data. 

B.   NOMENCLATURE 

The  following  parameters  and  arrays  were  used  in  the 
program;  their  respective  dimensions  are  shown  in 
parenthesis : 

ALV       array  used  to  store  interpolated  values  of 
dynamic  pressures  (feet) 

AMIXL     local  Prandtl  mixing  length  in  the  wake  (fee) 

CTAU      local  turbulent  shear  coefficient  for  each  point 
in  the  wake 

DELRAT    ratio  of  succeeding  DELSTR's  to  that  at  the  first 
data  station 

DELSTR    wake  displacement  thickness  (feet) 

DELX  axial  distance  between  succeeding  probe  tip 
positions  (feet) 

DIAMTR    diameter  of  wake  generator  (inches) 

EDDY      eddy  viscosity 

EPSIL  residual  in  momentum  integral  equation  when 
experimental  values  for  THETA  and  SHAPE  are 
substituted  into  it 

LABEL     label  for  subroutine  DRAW 

MN        number  of  data  points  laterally  through  the  wake 

MX        number  of  axial  data  positions 

MY        number  of  points  (always  odd)  into  which  the 
wake  is  divided 
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P         array  of  local  static  pressures  for  each  point 
in  the  wake  (psf) 

PATM  atmospheric  pressure  (inches  of  mercury) 

Q  freestream  dynamic  pressure  (psf) 

SHAPE  wake  shape  factor  (DELSTR/TKETA) 

THETA  wake  momentum  thickness  (feet) 

THETAP    inclination  of  static  pressure  manometer 
(degrees ) 

THETAV    inclination  of  dynamic  pressure  manometer 
(degrees) 

THRAT     ratio  of  succeeding  THETA' s  to  that  at  the  first 
data  station 

TITLE1    title  for  subroutine  DRAW 

TITLED    title  for  subroutine  DRAW 

U         array  of  velocities  for  each  point  in  the  wake 
(feet  per  second) 

URATIO    ratio  of  local  wake  velocity  to  local  freestream 
velocity 

XIN       array  of  probe  traverse  positions  through  the 
wake  (inches) 

XOUT      interpolated  probe  traverse  positions  for  each 
of  the  MY  points  in  the  wake  (inches) 

XPOS      probe  tip  position  measured  from  downstream 
side  of  the  wake  generator  (inches) 

YOUT      interpolated  manometer  readings  (for  YQ  and  YPS) 
for  each  of  the  MY  points  through  the  wake 
(inches) 

YPS       static  pressure  manometer  slant  distance  height 
difference  (inches) 

YQ        dynamic  pressure  manometer  slant  distance 
height  difference  (inches) 

MU(X)     kinematic  viscosity 

MM        index  of  the  maximum  velocity 
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MT 


index  of  the  maximum  velocity  in  the  upper 
half  of  the  data 


T  dummy  array  of  U  for  differential 

TRATIO  dummy  array  of  URATIO  for  differential 

ADP  dummy  array  of  ALP  for  differential 

D  dummy  array  of   for  calculations  Eddy  Viscosity 

NUM  array  of  MN 

MBM  array  of  MM 

MAD  array  of  MID 

MUT  array  of  MT 

MID  index  of  centerline 

C.   SUBROUTINES 


INTEG     uses  Simpson's  rule  to  numerically  integrate  the 
function  F.   It  must  be  supplied  with  K  points 
(always  odd)  equally  spaced,  DELX  units  apart 

QTFE   '   SSP3  library  subroutine  for  integration  using 
trapazoidal  rule 

QSF       SSP3  library  subroutine  for  integration  using 
Simpson's  Rule 

FILL      used  to  fill  in  arrays  to  maximum  dimension  with 
freestream  values 

PREIN     determines  midpoint  of  asignally  dimensioned 
array 

UTOP      determines  MT,  &  MM 

SPLIT     takes  the  entire  array  and  puts  the  values  above 
the  centerline  into  a  second  array  in  the 
input  parameters 


D.   USAGE 

The  details  on  the  use  of  this  program  are  listed  in  the 
comments  at  the  head  of  the  program.  The  deck  requires  120K 
to  compile  and  230K  to  run. 
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c  APPENDIX  E 

C 

£  DATA  REDUCTION  COMPUTER  PROGRAM 

Q  :Je  $  $  *!:  :je '!»  #:  S'  ^  *-?-'>>:  !    J-  -!'  :|<  *  >'<  #  ^  £  &  t<  %  ^  #  :$  it  $  £  £  4s  *  £  *  *  £  £  £  £  #  *  *  £  £  £  £  £  #  i^  * 

c 

C  A    HOT-WIRE    ANEMOMETER    STUDY    OF    FREE    TURBULENT    MIXING 

C  IN    AXIAL     PRESSURE    GRADIENTS. 

C 

c 

C  THIS    PROGRAM    WAS    DESIGNED    TO    REDUCE    P I TOT-STATI C , 

C  AND    HOT-WIRE    ANEMOMETER    DATA    TO    DETAILED    DATA 

C  CONCERNING    THE     PHYSICAL    PAPAMETERS    OF    THE    AIR    FLOW. 

C  IT    WAS     SPECIFICALLY    DESIGNED    FOR    USE     IN    THE    THESIS 

C  PROJECT    OF    LT.     M.     E.     DEARNEY     IN    1972,     HOWEVER,     THE 

C  PROGRAM    WILL    3E     USEFUL    TO    ANYONE     WHO    HAS    DATA    OF    A 

C  SIMILAR    NATURE. 

C 

C  TYPE    OF    DATA    REQUIRED 

C 

C  EITHER  PI  TOT  STATIC,  OR  HOT  WIRE  DATA  (WITH  A 

C  VELDCITY  REF.)  AT  EQUAL  INTERVALS  ACROSS  A  PARALLEL 

C  FLOW  2-OI MENS IONAL  WIND  TUNNEL.   THE  PROGRAM  IS 

C  PRESENTLY  SET  UP  TO  ANALYZE  3  SETS  OF  DATA  OF  12 

C  STATIONS  EACH,  THE  FIRST  BEING  PITOT  STATIC  VALUES 

C  AND  THE  SECOND  AND  THIRD,  HOT-WIRE  DATA.   ONLY  MINOR 

C  MODIFICATIONS  WOULD  BE  REQUIRED  TO  EXTEND  THIS  TO 

C  MANY  MORE  EXPERIMENTS. 

C 

C  INPUT: 

C 

C  1.  SIX  TITLE  DESCRIPTION  CARDS 

C  2.  ONE  LABEL,  INTERGERS  FROM  1-12 

C  3.  FIRST  RUN  TITLE 

C  4.  FIRST  2  LINES,  NUMBER  OF  STATIONS 

C  DIAMETER  OF  C YC LI NDE R { NOT  REQUIRED) 

C  5.  ONE  OR  MORE  CARDS  WITH  MANOMETER  ANGLES 

C  6.  ONE  OR  MORE  CARDS  WITH  ATMOSPHERIC 

C  PRESSURES. 

C  7.  TITLE   OF  FIRST  STATION,  FIRST  RUN 

C  8.  FIRST  3  LINES,  INTERGER  NUMBER  OF  POINTS 

C  K175. 

C  9.     REFERENCE    VALUE    OF    FIRST     POINT 

C  10.     DYNAMIC     PPESSURES    IN    I NC HE S ( PT-PS ) 

C  11.    STATIC    PRESSURES     IN     INCHES(PS) 

C  12.     SAME     DATA    7-11     FOR    FOLLOWING    STATIONS 

C  11.     SECOND    RUN    TITLE 

C  14.     NUMBER    OF    STATIONS 

C  DIAMETER    OF    CYCLINDER 

C  15.    SAME    AS    5,6,     £    7 

C  16.  NUMBER  OF  POINTS 

C  17.  REFERENCE  VALUE 

C  18.  VOLTAGE  DATA 

C  19.  THEN  CONTINUE  FOR  THE  REST  OF  THE  STATIC 

C  THE  THIRD  RUN  IS  THE  SAME  AS  THE  FIRST. 

C 

C  MAJOR    VARIABLES 

C 

C  U-y/ELCCITY 

C  CT4U-C0EFFICIENT    OF    STRESS 

C  MU(X),     VISCOUSITY    MODEL 

C  P-PRESSURE 

C  URATIC-    U/U(FREE    STREAM) 

C  ALD-    PSTATIC/12 

C  ALy/-    PDYNAMIC/12 

C  EDDY-    EDDY    VISCOUSITY 

C  AMISL-    REYNOLD'S    MIXING    LENGTH 

C 

C  MAJOR    SUBROUTINES 

C 

C  INTEG-    INTEGRATE    BY    SIMPSONS     RULE 


3.69 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

cc 

c 

c 

c 

c 


c 
c 
c 
c 


50 

51 

52 
53 
54 

55 
56 
57 
58 
70 
71 
72 
73 
74 
75 
76 
77 
78 
79 
80 


QTFE-  INTEGRATE  BY  TRAPAZOIDAL  RULE 

FILL-  FILL  THE  VELOCITY  AMD  PRESSURE  ARRAYS 

FOR  COMPUTATIONAL  PURPOSES. 
PRE'IN-  DETERMINE  THE  MIDPOINT  OF  THE  WAKE 

DEFECT. 
SM-  USED  IN  CONJUNCTION  WITH  LEAST  SQUARES 

ROUTINE  FOR  TAKING  DERIVITIVESt 

1=U,  2=DU. 
SPLIT-  DETERMINE  THE  MIDPOINT  OF  THE  FLOW 
UEE-  DETERMINE  THE  MAX  FLOW  VELOCITY 

A**£A£  £*&**:£$:  *  A  **•*£*  **H.'#£  £££**♦#  A  £  A*  #££&£****«£&**;!£&& 


IMPLIC 

REALMS 
REAL* 8 
REAL* 8 
REAL*8 
REAL*8 
REAL*8 
REAL* 4 
D I M  E  N  S 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMENS 
DIMFNS 
COMMON 
DATA  G 
H=0.1E 
Y=H/12 


I  T  R 
TI 
TI 
TI 

LA 
LA 
LA 
LA 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

ION 

/A/ 

A  MA/ 

0 

.0 


EAL- 
TLE1 
TLE2 
TLE4 
BEL1 
8EL2 
BEL3 
BEL( 
YFAT 
ALP( 
AMU( 
A(  12 
EDO/ 
URAT 
UR(1 
YOUT 
AXIS 
DELS 
T(l 
0(  12 
THET 
WX(1 
XFIT 
RAKE 
FAKE 
NUM( 
H 
62.3 


4     (A 

(  12) 

(12) 

(  12) 

/' 

/« 

/  ' 

12) 

(12) 

12,1 

12)  , 

,  175 

(12, 

I  3(1 

75), 

(175 

(  175 

TR(1 

2,17 

)  ,  DE 

AP(1 

2 )  ,  W 

(  12) 

(175 

(  175 

12) 


-H,0-Z) 


DELS'/ 
THET'/ 
SHAPV 

, YFIT( 12) ,YFOT( 12) ,YFUT(12) 
75),ALV(12,17  5),U(12,17  5),P(12,175) 
MUT(12) 
),B(12,175),C(12,175),D(12,175) 


175)  , AM  I  XL (12, 17  5),X0UT(  12, 
2,175) ,CTAU( 12,175) ,YQ( 175) 


UD(175) 


175) 

,YPS(175) 
,F(175) 


UP(  175) ,UC(  175) 

)  ,XIN(175) 

) 

2),THETA(  12),  SHAPE  ("12)  ,  EPSILU2) 

5)  ,ADP(12,175) 

LRAT  (  12)  ,  THRAK  12)  ,THETAV(  12) 

2 ),  PATM(12) ,XPOS( 12),DELY{ 12) 

H(12) , IW(12) 

) 

) , MAD (12) ,MBM(12) 


2/, RHO/0. 002378/ 


SET     INPUT/OUT>UT    FORMATS 


FORM 
1'  ) 

FORM 
FORM 
FORM 
FORM 

l1) 
FORM 
FORM 
FORM 
FORM 
FORM 
FORM 
FORM 
FORM 
FOR^i 
FORM 
FORM 
FORM 
FORM 
FOR^ 
FORM 

l'STA 


AT(  • 

AT(  I2,F10.2) 
AK8F10.1) 

AT (8F10.2) 

AT(  ' 

AT( 13) 

ATOF10.3) 

AT(///////«  ' 

AT(///////'  ' 

AT(////////// 

AT( //» 

AT (//' 

AT( //' 

AT(//« 

ATI/// 

AT(//« 

AT(//« 

AT(//' 

AT(  '  1' 

AT  (•  1' 

TIONS' 


,T20, 
,T20, 
,  T  2  0  , 
,T20, 
'  ,T53 
,T10, 
,T10, 
T10 


T20,6A8) 
T20, 12A4) 
•  -T33  ,«  ECHO-CHECK  DATA'///) 
('  •  ,T33,8F10.2//) ) 
(•  ',T33, 6F10.2//) ) 
'  ,T33,8F10.2//)  ) 
• ,T33,8F10.2//) ) 


THETAV 
THETAP' 


PATM 
XPOS 
«MN  = 
XIN« 
YO  • 
YPS' 


( 

13///) 

•  '  ,T21 ,8F10.3//)  ) 

) 

) 


,T21,8F10.3//) 
,T21,8F10.3// ) 


,T50, '**-RFSULTS  OF  DATA  REDUCTION***') 
//////'  ',T50,  'AXIAL  STATION'  , I  2,  'OUT  OF'  , 12 
/////) 
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HRMAK •     • ,T18, 'DEL STAR=« ,G 12.5, 5X, «THETA=» ,G12.5,5X, 
JcuapE=' ,G12.5,5X, »EPSILON=«,G12.5) 
]:ORMAT(  //////'     ?*T8,«  POINT'  ,T30  , «  U'  ,T48  ,  •  U/UE"  ,  T70 ,  •  P« 

1:0RMAT^0^  !t9  ,  13  ,T24,G1  2.  5,  T44 , G12 .5 ,T64, G12. 5,T  84, 

loRMAT( • l',T44, '***SUMMARY    OF     INTEGRATED    WAKE    PROPER', 

toRMATC/////'     '. T3, 'AXIAL*  ,T75 , • THETA' / "     •  ,T2 • •■ ST ATION 
•  <  .T16     'X/D«     T2  9,  'DELSTAR«,T44,  • DELRAT 10 ' ,T60, ' THETA '  , 
T75,  'RATIO'  ,T90,  'SHAPE'  ,  T105 , ' E PS  I LON ' )  _    Tirz 


=0RMAT(//«     ' ,'  LAST=«  ,13) 

[BgSSiu*>b?>/?A"Te,;??oriffjT3§'.,U',T48..U/UE..T70, 

,T88, 'GTAU'  ) 

;-8rMAT('0'(t60,« TURBULENCE    CHE CKS' // /Tl 0  ,'  STAT  ION « , 

15X. 'ETA', 12X, 'MU(X)'//)  „     #/1 

r  FORMAT  t'     '  ,  T14  ,  I  2  , 1 IX  ,  Fl  2  .4  ,  5X  ,  Fl  2  .4/  ) 

READ(5,90)  TITLE1 

READ(5,90)  TI TLE2 

READ(5,90)  TITLE4 

READ(5,92J  LABEL 

ZERO    ALL    ARRAVS     BEFORE    COMPUTATION    BEGINS 

THIS    LOOP    CARRIES    THE     PROGRAM    THRU    THREE    SETS    OF    DATA 


**#********#****************' 


DO    1000    KK=1,3 

Z  =  D.D 

DO    1    1=1,12 

DELRAT( I )=Z 

DELSTRU  )=Z 

DELYU  )=Z 

EPSILt I )=Z 

PATM1I ) =Z 

0(1  )  =  Z 

SHAPE ( I )=Z 

THETAl*  I  )=Z 

THETAP1 I )=Z 

THETAVt I )=Z 

THRATC I j=Z 

AMU( I )=Z 

XPOSCI  )=Z 

DO    1    J=l,175 

A(I, J)=Z 

B(T , J) =Z 

C(I ,J)=Z 

CTAU( I, J)=Z 

D(I , J)=Z 

F( J)=Z 

P( I, J)=Z 

URATIOC I, J )=Z 

X0UT( I , J)=Z 

GO    TO     (36, 1, 1) ,     KK 

ALP(I , J)=Z 

ALVII , J) =Z 

U( I, J)=Z 

CONTINUE 
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J7=MZ 
MZ=MZ+1 

DO    864    J=MU,17  5 
YPSU)  =  ALP(I,ME)*12. 

YPSU7)=YPS(MZ  ) 

IF( J7.EQ.1)     GO    TO    862 

J7=J7-1 
8  64  CONTINUE 
862  CONTINUE 

GO  TO  4  3 
41  READ{5,56)    ( Y PS { J )  ,  J  =  l , MN ) 
C 

C      ECHO-CHECK  INPUT 
C 

43  WRITE(6,95) 

WRITE (6, 54) 

WRITE(6,75)   MN 

WRITE (6* 76)   ( XI Nf J)  ,J=1,MN) 

WRITE{6,95) 

WRITE{6,77)   (YQ(J)  ,J  =  1,MN) 

WRITE{ 6,95) 

WRITE (6, 78)   (YPS( J  ),J=1,MN) 

MY=MN 

NUM( I )=MN 
C 

C      CALCULATE    U  VELOCITY 
C 

IFtKK.GT.l )  G3  TO  37 

DO  3  J=1,MN 

YOUT( J)=YO(J  ) 

ALV( I, J)=YOUT( J )/12.0 

U( I , J)= SORT (2. 0*GAMA*SIN{THETAV(I) }*ALV( I , JJ/RHO) 
3     CONTINUE 

GO  TO  38 

37  MM=M3M{ I  ) 
UDUM  =  U(  I, MM) 

CALL  PREINtYOt MM, MID) 

MM=MID 

DO  632  J=1,MY 
632  U(I , J) =YO( J) 

UTOP=UEE<  U,M^,UMAX,MM,I ,MT) 

MV=MBMl  I  ) 

ADUM=ALV(I ,MV) 

DO  39  J  =  1,MY 

ALV(I,J)=(      ADUM*YQ(J))/YO(MM) 

U{ I , J) =( UDUM*YQ( J) ) /YO(MM) 
39  CONTINUE 
C 

C      DETERMINE  MAXIMUM  FREESTREAM  VELOCITY 
C 

c 

38  CALL  PREIN(YO,MN,MID) 
MM  =  MID 

UTOP=UEE  (U,MN,UMAX,  MM,I,MT) 

0( I )=0. 5*RHO*J( I ,MM)**2 

WRITE(6,55)  MID 

WRITE (6,  55)  MM 

WRITE{6,55)  MT 

WRITE(6,53)  UMAX 

WRITE{6,53)     UTOP 

MAD( I)=MID 

MBM(I)=MM 

MUT( I ) =MT 
C 

C      CONTINUITY  CHECK 
C 

DO  99  J=1,MY 
99  F( J)=U( I ,J) 

J4=2 

CALL  SPLIT  (F, FAKE, MT, MID,  J4) 

CALL  INTEGtFAKE, J4,Y, RESULT) 

WRITE(6,93)  RESULT 
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c 

C      CALCULATE  DISPLACEMENT  THICKNESS 
C 

DO  3  1  J  =  1,MY 

URATIOC  I  ,  J)=Ul  I, J  J /U( I, MM) 

F( J) =1. 0-URATI 0(1 , J) 
31  CONTINUE 

Jl=2 

CALL    SPLIT(F,FAKE,MY,MID,J1) 

CALL  QTFEfY, FAKE, RAKE, Jl) 

DELSTRt  I  )  =  RAKEU1  J 

DELRATd  )=DELSTR(I  )  /DELSTR(l) 
C 

C      CALCULATE  MOMENTUM  THICKNESS 
C 

DO  4  J=1,MY 

F( J)=URATIO( I, J )*(1.0-URATIO(  I, J)  ) 

4  CONTINUE 
J2  =  2 

CALL  SPLIT (F, FAKE, MY, MID, J2) 

CALL  OTFE(Y, FAKE, RAKE, J2) 

THETAC I >=RAKE( J2) 

THRATt I J=THETA(I )/THETA(l) 
C 

C      CALCULATE  SHAPE  FACTOR  IN  MOMENTUM  INTEGRAL 
C 

SHAPE  (  I  )  =DELSTR(  I  )  /THETAU  ) 
C 

C      CALCULATE  THE  PRESSURE  FIELD 
C 

DO  5  J=1,MY 

IF(KK.GT.l )  GO  TO  40 

X=XOUT( I  ,J) 

YOUT( J)=YPS( J) 

ALPl I, J)=YOUT( J J/12,0 

F( J)=URATTO{ I , J) 
40  P( I, J)=PATM( I)-  GAMA*SIN(THETAP(I) )*ALP( I , JJ 
C 
C 

5  CONTINUE 
C 

c 

CALL  DRAW(MN, XIN,F, 0,0, LABEL, TITLE1, 0,0,1  ,1,2,2,6,  6,1 
C 

c 

6  CONTINUE 
C 

c 

C  FILL     IN    FUNCTION    ARRAYS    TO    MAX    DIMENSION    FOR 

C  CALCULATIONS 

C 

c 

CALL    FILL (NUM. 1,U) 

CALL    FILL! NUM, 1 ,ALP)     , 
C 

C  PREPARE    DUMMY    VARIABLES    CENTERED    ON    THE    FLOW 

C  CENTERLINE 

C 

DO    711    LTB=1,2 

IF(LTB.EQ.2)     30    TO    712 

DO    451    1=1, MX 

MID=MAD( I) 

DO    451    J=MID,175 

JK=J-MID+1 

T( I  ,JK)=U( I , J) 

ADP( I , JK)=ALP( I , J) 
451    CONTINUE 

GO    TO    714 
712    DO    713    1=1, MX 

MID=MAD( I) 

JL=MID    +1 

DO    713    J=MID,100 
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c 

c 
c 


713 
714 


C 

c 

c 
c 
c 
c 


c 
c 

c 


575 
574 

573 

576 


JK=J-MID+1 

JL=JL-1 

T( I t JK)=U( I , JL ) 

ADP(  I, JK)=ALPi  It JL ) 

IF (JL.EO.l)     GO    TO    714 

CONTINUE 
CONT INUE 

CALCJLATE    RESIDUAL    IN    MOMENTUM    INTEGRAL    EQUATION 


1=1 

DELX=X 

MV=MBM 

MM=M3M 

E  1  =  (  T  H 

E2=(TH 

E3=(U( 

EPSILt 

DO    7     I 

DELX=X 

M1=M3M 

M2=MBM 

M3=MBM 

MY=NUM 

E1=(TH 

E2=THE 

E3  =  { U  ( 

EPSIL( 

CONT  IN 

I  =  MX 

DELX=X 

MV=MBM 

MM=MBM 

E1=(TH 

E2=(TH 

E3=(U( 

EPSIL( 


POSt 
(  1+1 
(I) 
ETA( 
ETA( 
1+1, 
I  )  =  ( 
=  2  , M 
POS( 
(  1-1 
(  I  ) 
(  1  +  1 
(  I  ) 
ETA( 
TA(  I 
I+lt 
I  )=( 
UE 


i+n-xposu  j 

) 

I+1)-THETA(I ))/DELX 

1  )*( SHAPE( I  )+2.)  ) 

MVJ-U(  I  f  MM)")  /(DELX*U(I,MM)  ) 

E1+E2*E3)/(E1-E2*E3) 

X  1 

H-1J  -XPOS(I-l) 

) 

) 

1+1) -THETA(I-l) )/DELX 

)*(SHAPE( I )+2.0) 

M3)-U( I-1,M1 ))/(DELX*U( I,M2) ) 

E1+E2*E3)/(E1-E2*E3] 


POS( I )-XPOS(I-l) 

{ I ) 

( 1-1 } 

ETAU  )-THETA(  1-1  )  )/DELX 

ETA( I )*< SHAPE (I )+2e) ) 

I,MV)-U(  I-1,MM))/(DELX*U( I,MV)  ) 

I }=IE1+E2*E3)/(E1-E2*E3) 


CALCULATE    LOCAL    TURBULENT    SHEAR    COEFFICIENT    IN    WAKE 
STATION    ONE    -    FORWARD    DIFFERENCE 


1=1 

MM  =  MBM(  I  ) 

MY  =  NUM(  I )-MAD( I )-K 

DELX  =  XPOS(  I  +  l)-XPOS(I  ) 

DO    576    K=l , MY 

DO    575    J  =  1,K 

A(  I, J)  =  2.*T(I, J)-T(  I,K J 

B(I t J) =(T(I+1, J)-T( I* J) )/DELX 

C(  I,  J)  =  (ADP(IH,  J)-ADP{I  ,J)')/(DELX*ALV{  I  ,MM)  ) 

CONT  INUE 

J4  =  0 

DO     574    J  =  1,K 

J4=J4+1 

F( J)=A(I  ,J)*B(  I , J) 

CALL    QTFE( Y,F, RAKE, J4) 

J5  =  0 

DO    573    J=1,K 

J5=J5+1 

F( J)=C( I, J  ) 

CALL    QTFE( Y,F,FAKE, J5) 

CTAUt I,K)=(2./U( I,MM)**2)*RAKE( J4)-FAKE( J5) 

CONTINUE 

STATION!    2    TO    (MX-1)     CENTRAL    DIFFERENCE 

DO    10    1=2, MX1 

DELX  =  XPOS( I  +  1)-XP0S(I  ) 

MM=MBM( I ) 

MY=NUM( I )-MAD( I )+l 
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c 
c 
c 


M I  D=  1 

DO  10  K=MI 
DO  9  J=MID 
MA=J-MID+1 
A( I , J)=2.* 
BCIt J)-(T( 
C(  I,  J)  =  { AD 
9  CONTINUE 
J  4=0 

DO  701  J=M 
MA=J-MI  D+l 
J4=J4+1 

7  01  F(MA)=A(I, 
CALL  OTFE( 
J5=0 

DO  70  2  J=M 
MA=J-MI D+l 
J5=J5+1 

7  02  F(MA)=C(I* 
CALL  QTFF! 
CTAU( I ,K)= 
10  CONTINUE 
IF{ LTB.EO. 
DO  521  1=1 
MY=NUM( I )- 
MID=MAD(  I  ) 
DO  52  1  J=l 
JK=J+MID-1 
D( I ,JK)=CT 

5  21     CONTINUE 
GO    TO    711 

718  DO  719  1=1 
MY=NUM! I }- 
DO    921    J=l 

921    F!J)=CTAU( 
MY=NUM( I ) 
L=MAD! I 
DO    920 
L  =  L-1 
CTAU( I , 
IFCL.EO. 1) 
CO NT INUE 
L=MAD( I  ) 
DO    923    J=L 
CT  AU  !  I  ,  J  )  = 
CONTINUE 
CONTINUE 
WRITE(6,79 
DO    14    1=1, 
WRITE (6, 80 
WRITE(6,81 
WRITE!  6,  82 
MY=NUM( I ) 
DO    13     J=l 
IF( J.EO 
IF! J.FQ 
IF(  J.EO 
IF( J.EO 
WRITEC6 

13  CONTINUE 

14  CONTINUE 
WRITE!6,84 
WRITE(6,85 
DO  15  1  =  1, 
XD=12.0*XP 
WRITE! 6,86 

1 , SHAPE!  I  ), 

15  CONTINUE 


D ,  MY 
,K 

T(I , J)-T  !I,K) 

I+lf J)-T(I-1,J) )/(2.*DELX) 

P( 1+1, J  J-ADP!  1-1, J ) )/(2.*DELX*ALV( I, MM)  ) 


ID,< 


J)*B{ I  ,J  ) 
Y,F, RAKE, J4) 

ID,< 


J) 

Y,F,FAKE,J5) 

( 2./U! I,MM)**2)*RAKEI J4)-FAKE( J5J 

2)  3  0  TO  718 
,  MX 

MAD! I)+l 

,MY 

AU!I ,J) 


MX 
MAD! I )+l 
,  MY 
I  ,J) 


) 

J=2,100 


L) 


920 
922 

923 
719 
711 


30 
60 
90 
12 

83 


F(J) 

GO  TO  922 


,MY 

D!  I,  J) 


) 

MX 

)   I  ,  M  X 

)  DELSTR!  I  },THETA(  I  },  SHAPE!  I),  EPSILU) 

) 

MY 

)  WRITE  I  6, 94) 

)  WRITE(6,94) 

)  WRITE (6, 94) 

0)  WRITE!6,94) 

)  J,U(I,J),URATIO!I,J),P!I,J),CTAU(I,J) 


) 

) 

MX 

OS! I J/DI AMTR 

)  I  ,XD, DELSTR (I )  ,DELRAT(I)  ,THETA(I) ,THRATII) 

EPSIL  (  I  ) 


CALCJLATE  THE  TURBULENT  KINETIC  ENERGY  RESIDUAL  AND  MU 
DO  715  1=1, MX 
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715 


231 
230 


2  36 


233 


232 


237 


581 


MY=NUM( 

MID=MAD 

DO    715 

JK=J+MI 

CTAU( I  * 

CONT1NU 

DO    230 

J4  =  0 

MID=MAD 

MT=MUT t 

DO    231 

J4=J4+1 

MA=J-MI 

F  (  MA )  =  ( 

CALL    OT 

YFIT( I ) 

CONTINU 

MT=MUT( 

MID  =  1 

J4  =  0 

1=1 

DELX=XP 

MY=NUM( 

MID=MAD 

DO    236 

MA=J-MT 

F(MA)=C 

J4=J4+1 

CALL    OT 

YFOTI 

YFATI 

YFOT(  I  ) 

DO    232 

MT=MUT( 

DELX=XP 

J4=0 

MID=MAD 

DO    233 

M.A  =  J-MI 

Ft MA)=C 

J4=J4+1 

CALL    OT 

YFOTI 

YFATI 

YFOTt I  ) 

CONTINU 

MT=MUT( 

J4  =  0 

I=MX 

MID=MAD 

DO    237 

MA=J-MI 

F(MA)=C 

J4=J4+1 

CALL     OT 

YFOTI 

YFATI 

YFOTt  I  ) 

1=1 

MT=MUTt 

DELX=XP 

MID=MAD 

J4  =  0 

DO    581 

MA=J-MI 

FtMA)=U 

J4=J4+1 

CALL    OT 

DER=(YF 

AMUt  I)  = 

DO    234 

DELX=XP 


I )-MAD( 

tl) 

J=1,MY 

D-l 

J)=D(I, 

1=1, MX 

{  I) 

I  ) 

J  =  MID,M" 

D+l 

U(I ,MT) 
FEt Y,F, 
=  P,AKE(J 
E 
1) 


I  )  +  l 


JK) 


*=*2-U(I  , J)**2)*U(I, J) 

RAKE , J4) 

4)~RH0 


OSt  1+1) -XPOS  tl  ) 

I  ) 

t  I) 

J=MID,MY 

D+l 

TAUII,MA)*2.*((U(I,J-1)-U(I,J+1))/(2.*Y)) 


FEt Y,F, 

=RAKE(J 

=RAKE( J 

=YFOTI/ 

I=2,MX1 

I  ) 

OSt 1+1) 

(  I) 
J=MID,M 

D  +  l 
TAUt  I  ,H 

FE(Y,F, 

=RAKE( J 

=RAKE(J 

=YFOTI/ 

E 

MX) 


^AKE, J4) 

4)-(YFIT(I)-YFIT(I+l))/DELX 
4)+t YFlTtI J-YFIT t 1+1) )/DELX 
YFATI 


XPOStI  ) 


Y 

A)*2.*l (U(I  , J-l)-U(IfJ+l) )/ (2.*Y) J 

RAKE, J4) 

4)-(YFITt 1-1 J-YFIT t 1+1) )/(2.*DELX) 

4)+( YFITt I-l)-YFITt 1+ 1 ) ) / ( 2 .*DE LX) 

YFATI 


(  I) 
J=MIDtMY 

D+l 

TAUt  I,MA)*2.*(  (Ut  I, J-l)-U(  I, J+l)  ) /t2.*YJ) 

FEt  Y,p, RAKE, J4) 

=  RAKE(J4 )-(YFIT  (I-l)-YFIT(I  ))/DELX 

=RAKE(J4)+( YFIT (I-l)-YFIT t I) )/DELX 

=YFOTI/YFATI 

I  ) 

ost  i  +  ij-xposu  j 
tn 


J=MID,M 

D  +  l 

( I, J)** 

FE(Y,F, 
I TC I+1J 

DER/RA< 
I  =2  ,MX1 
OSt 1+1) 


RAKE, J4) 
-YFITtI ) J/DELX 
E(J4) 

-XPOStI ) 
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MT=MUT(  I  ) 

MID  =  MAD(  I  J 

J4=0 

DO    235     J=MID,MT 

MA = J- MI  0  +  1 

F(MA)=U( I, J)** 2 
235     J4=J4+1 

CALL    OTFE(  Y,F,RAKE,J4) 

DER=(YFIT( I+1)-YFIT{I-1))/(2.*DELX) 

AMUU ) =  DER/RAKEU4) 
234    CONTINUE 

WRITE(6,96 ) 

DO    717     1=1, MX 
717    WRITE(6,97)     I , YFOT ( I )  , AMU(  I  ) 
C 

C  CALCULATE     EDDY    VISCOSITY    AND    REYNOLD'S    MIXING    LENGTH 

C 

WRITE (6, 87 ) 

DO    17    1=2, MX 

WRITE(6,88)     I 

MIN=(  (NUM(  I )-l  )/2)  +  10 

MAX=NUM( I )-10 

DO    16    J=MI N.MAX 

E1=ABS(  DU,J)*Q(I)) 

E2  =  ABS(U(I  ,  J  +  1J-UU  vJ-1) ) 

E3=2.0*H 

EODY(  T  ,  J)  =  (E1/  RH0)/(f:2/E3) 

AMI  XL  (I  ,  J)  =  SQ^T(  E1/(RH0*(E2/E3)**2)  ) 

WRITE (6, 89 )       J , EDDY(  I, J  )  ,AMIXL( I, J) 

16  CONTINUE 

17  CONTINUE 
C 

1000    CONTINUE 
999    STOP 
END 
C 
C 
C 
C 
C 

Q  !>s  f,;  £  j£  #  £  #  £  iSc  *■  *•  %:  >fc  sfc  ^cs(<  fc  £  £  >)c  :!t  :£  #  sfc  £  #  &  $  *!f  £  ^  &  £  &  #  ^  #  ^  £  *  ^  ^  #  #  ^^  s? '?'  £ 

c 
c 

C  THIS    PROGRAM    WAS    DESIGNED    TO    INTEGRATE    BY    THE    USE    OF 

C  SIMPSONS    RULE. 

C 

C  INPUT: 

C 

C  CALL    INTEG(F,M, DELX, RESULT) 

C 

C  F-  ARRAY  OF  VALUES  OF  THEINTEGRAND 

C  M-  THE  MUMBER  OF  POINTS 

C  DELX-    THE    STEP    SIZE 

C  RESULT-  THE  RESULT  OF  THE  INTEGRATION 

C 

c 

c 

c 

SUBROUTINE     INT  EG ( F , K, DELX, RESULT ) 
IMPLICIT    REAL*4    (A-H,0-Z) 
DIMENSION    F(l) 
C  INTEGRATION    BY     SIMPSONS    RULE 

Kl-K-1 
A=F(1)+F(K) 
B=D.O 

DO     1     I  =  2  ,  K  1 ,  2 
SUM=4.0*F( I ) 
B=B+SUM 
1     CONTINUE 


178 


C=0.0 

DO  2  I=3,K,2 

SUM=2.0*F< I) 

C=C+SUM 
2  CONTINUE 

RESULT=DELX*U  +  B+C)/3.0 

RETURN 

END 
C 
C 

c 
c 

C         THIS  SUBROJTIME  WAS  DESIGNED  TO  FILL  AND  ARRAY  WITH 

C         VALUE  EXISTING  AT  AN  INDEX  IN  THE  ARRAY. 

C 

C         INPUT:      CALL  FILL(M,N,U) 

C 

C  M-    AN    ARRAY    OF    INTERGER    VALUES    OF    THE     INDEX 

C  N-    A    FUNCTION     VARIABLE 

C  U-    A    DOUBLE    DIMENSIONED    ARRAY    OF    VALUES,     U(I,J) 

C  WHERE    THE     INDEX    M    IS    CONTAINED    IN    J. 

C 

C 

c 
c 

SUBROUTINE  F I LL ( NUM , MX » U) 

IMPLICIT  REALM  (A-H,0-Z) 

DIMENSION  NUM( 1) ,U(12,1) 

DO  2  J=1,MX 

M=NUM( J) 

DO  1  I=M,175 

1  U( J.  I)=U(J,M) 

2  CONTINUE 
RETURN 
END 

C 

c 
c 

c 
c 
c 

C         THIS  SUBROUTINE  WAS  DESIGNED  TO  TAKE  AN  ARRAY  AND 

C      DETERMINE  ITS  MID-POINT  BY  MEANS  OF  AN  ARBITRAY 

C      FUNCTION,  AND  RETURN  THE  INDEX  OF  THAT  POINT. 

C 

C  INPUT: 

C 

C  CALL  PREIN(A,N,M) 

C 

C  A-  THE  ARRAY  OF  VARIABLES 

C  N-  THE  LIMIT  OF  THE  VALUES  TO  BE  CHECKED 

C  M-  THE  INDEX  OF  THE  MID-POINT 

C 

C 

C 
C 

SUBROUTINE  PREIN(A,M,M) 

DIMENSION  A(l) 
C 

C  DETERMINE    MIDPOINT    OF    MAKE    DEFECT 

C 

DUM=1000. 

Nl=N-2 

DO    1    1=3, Nl 
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B=A( I-2)+A(I-l )+A( I )+A( 1+1 }+A< 1+2) 
C=AMIN1(B,DUM) 

IF(C.GE.DUM)    GO    TO     1 

DUM=B 

M=I 

1  CONTINUE 
RETURN 
END 

C 
C 

c 
c 

C        THIS  SUBROUTINE  WAS  DESIGNED  TO  TAKE  AN  ARRAY  OF 

C      VALUES  AND  SPLIT  THEM  AT  AN  ARBITRARY  POINT,  AND 

C      PLACE  THESE  VALUES  IN  ANOTHER  ARRAY. 

C 

C  INPUT: 

C 

C  CALL  SPLITC A»B»M,Nt JK) 

C 

C 

C  A-  THE  INPUT  ARRAY 

C  B-  THE  OUTPUT  ARRAY 

C  M-  THE  LIMIT  OF  THE  SPLIT 

C  N-  THE  START  OF  THE  SPLIT 

C  JK-  THE  NUMBER  OF  VALUES  BETWEEN  M,  AND  N 

C 

C 

c 

c 

SUBROUTINE   SP L IT ( A  ,B  ,  M , N, JK ) 
DIMENSION   A(I ) ,P( 1) 
IF(JK.CQ.l)  G3  TO  1 
DO  2  I=N,M 
B(I-N+1)=A( I ) 

2  CONTINUE 


JK=M-N+1 

GO  TO  5 

1 

MD=M-N-1 

DO  3  I =1,MD 

IJ  =  N-I4  1 

B(  I)=A(  IJ) 

IF(I J.LE.l) 

3 

CONTINUE 

4 

JK  =  MD 

5 

RETURN 

END 

GO  TO  4 
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c 
c 
c 

c 
c 
c 

C        THIS  FUNCTION  WAS  DESIGNED  TO  SELECT  THE  LARGEST  VAL 

C      ARRAY*  AND  THE  LARGEST  VALUE  IN  THE  TOP  HALF  OF  THAT  A 

C 

C 

C  INPUT: 

C 

C 

C  FUNCTION  UEE(U,MY,UM,MM, I,MT) 

C 

C  U-  THE  ARRAY  TO  BE  SEARCHED 

C  MY-  THE  EXTENT  OF  THE  ARRAY  TO  BE 

C  SEARCHED,  AN  INDEX 

C  UM-  THE  MAXIMUM  VALUE  IN  THE  ARRAY 

C  MM-  THE  INDEX  OF  THAT  VALUE 

C  I-  THE  INDEX  OF  THE  FIRST  DIMENSION 

C         ,  0  F  U  (  I  ,  J  ) 

C  MT-  THE  INDEX  OF  THE  LARGEST  VALUE  IN 

C  THE  UPPER  HALF  OF  THE  ARRAY 

C  UEE-  THE  VALUE  OF  THE  LARGEST  VALUE  IN 

C  THE  TOP  HALF  OF  THE  ARRAY 

C 

c 

c 
c 

FUNCTION  UEE(J,MY,UM, MM, I , MT J 

DIMENSION  U(12,l) 

UM  =  0. 

M=MM 

DO  1  J=1,MY 

UI=UM-U( I, J) 

IFCUI )  2,2,1 

2  UM=U( I , J ) 
KM  =  J 

1  CONTINUE 
UEE=0. 

IF(MM.GT.M)  GO  TO  3 
DO  4  J=M,MY 
UI=UEE-U(I  ,J  ) 
IF(UI)  5,5,4 

5  UEE=U( I , J) 
MT=J 

4  CONTINUE 
GO  TO  6 

3  UEE=UM 
MT=MM 

6  RETURN 
END 


THE    FOLLOWING     IS    A    PRINT    OUT    OF    THE    RAW    PITOT-STATIC 
DATA    FOR    THE    STUDY 


//G0.FT06F001     DD       S P ACE= ( CYL ,  ( 1 ,  1  )  ) 
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